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ON THE ANALYTICAL FORM OF THE EARTH'S MAGNETIC ATTRACTION 
EXPRESSED AS A FUNCTION OF TIME 

V. Carlheim-Gyllenskold** 

A complete theory of the Earth's magnetic attraction /3* 
must include the force's manifestations for. 'any given time. 

The old theory solves this problem for a specific per- 
iod of time, but does not account for an aspect called 
secular variations. Generally, the variations of elements 
observed in a given place are expressed by empirical for- 
mulas, using either a series operating according to the 
powers of time, or a trigonometric series. However, this 
does not teach us anything about the nature of things, since 
it is a well known fact that any analytical function is 
expandable, within a finite domain, into a Taylor series 
or into a Fourier series. Furthermore, these methods are 
totally useless for preliminary calculations. 

As the variations are highly complex in appearance, 

we naturally decide to look for simple laws, if there are 4 

! 

any, not in the effect observed, but in the cause of the I 

phenomena, i. e. in the distribution of the magnetic mass ,,j 

inside the earth. This distribution is expressed by arbi- d 

\) 

traries shown in the expansion of the function of the i 

earth's magnetic forces which precisely are the integrals 
defined in certain functions of the mass coordinates. It 
is clear that any theory must commence by a study of these 
arbitraries. 

This problem, seemingly very simple, in reality en~ 

t 

counters serious problems due to the nature of the question 
itself, but also to the imperfection of the old observations. 

If the illustrious author himself of the Allgemeine Theorie 

des Erdmagetismus despaired in being able to succeed with J r . 

v 

*Numbers in the margin indicate pagination in the original text. 

^Stockholm 1896 
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with the materials bequeathed to us by our predecessors } 
undoubtedly, it was a sort of scientific gourmandism, 
not too commonly foijnd among scholars. Bessel is even 
more skeptical about this in his PopulUre Vorlesungen 
Uber wissenschaf tliche GegenstMnde (P.339); 372). All 
contemporary authors have respected the authority of 
these scholars, and for good reason. Spoiled by the ever 
increasing precision of modern observations, we are too 
inclined to look down on traditional observations made 
with imperfect instruments, by unskilled observers, yet 
when they are combined in the appropriate manner, they 
are often called upon to render their good services. 

The following lines intend to show that this task 
is possible, by presenting a simple method by which the 
objective can be reached without too much difficulty, 
and by presenting the already very satisfactory results 
that I obtained in the first application of this method. 

As we were confronted only with unknowns in the sub- 
ject in question, we had to begin by looking for the em- 
pirical laws of the phenomena by inductive reasoning. 

This is the first aspect that we will deal with. 

Prom this, by elementary reasoning, we will deduce the /4 
magnetic field which produces each instant change in the 
earth's magnetization, then we will trace back to the phys- 
ical causes which we think are responsible for the results 
obtained. Once we have found this principle, we will be 
able to express the differential equation that the solution 
of the question depends upon and we will see the main 
features revealed by experience follow from this at once. 
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ORICMNAL PASS IS 

OF POOR QUALITY 


The function of the forces of a sphere of attrac- 
ting material, of radius R, for an outside point of dis- 
tance v from the center of the sphere, may be expanded 
in a decreasing series with respect to the powers of v 
and therefore having this expression: 


V 


_ 1 Jo) 

~ v' 1 



.Y (1 } 



+ etc . 


where the Laplace functions are represented. We know that 
Y ^ is a rational and whole function of / u, rr, 

and VT — .costs. where /« = cos 0, which the vector radius r 
forms with the x, ra 1 . axis. The is the angle which the plane 
passing through the radius and through this axis forms with 
the x and y plane. 


(i J 

Each of the Y functions contains 2i+l arbitrary 
constants, each term is made up of 2i+T terms. The part 
of Y ^ depending on the »ts angle will be 

I^ 0 = (i — y-y> . . (A™ sin net + JB ^ 0 cos nis). 


by letting stand for the Legendre* polynomials. If 

we do n = 0 3 n = 1, n = 2 3 ... n = i 3 in this function suc- 

sessively, the sum of all the resulting functions will be 

the general expression of Y^K The A ^ and B^ are the 

n n 

numerical coefficients which are dependent upon the mass 
distribution and which must be derived in the experiment. 


We are proposing an empirical study of the secular 

inequalities of arbitraries A {% and B , based on acces- 

n n 

sible data. 


Following the usual procedure, we will divide the 


*We have retained the notations of la Mecanique Celeste 

[Celestial Mechanics] (book III, chapter .II); if we wanted . 

to use Gauss' notations, we would set A'' 1 ' ~ h 1 ' ,n , = g' l,n 

n n a 
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the attracting force S, for any given location on the 
earth's surface, into three rectangular axes, by adopting 
the plane of: the horizon for the xy plane and the vertical 
plane for the z axis. Once this is established, the force 
components will be the partial derivatives of V with res- 
pect to these axes: 


v ,/T 5 dV 


Y. 


dV 


y i—pi'da' 


dV 

W' 


The components of the attraction are linear functions of 

/ * 1 / • » 

the arbitraries A^ ' and which may be determined if 

a sufficient number of data are used, for any time where 
the force is known. 


The 'criteria were different for past centuries. 

Earlier observations did not give us the magnitude of the 
magnetic force, but only its direction. It is therefore /5 

necessary to establish equations dependent upon the declin- 
ation and the dip of the magnetic needle and to select an 
expression from Which the unknown values may be deduced. 


The most direct method of obtaining the differences 

(i) (i) 

from the -.constants A and B with respect to time 

n n e 

would be to calculate them using the formulas which give 
the differences of the declination and the dip as a func- 
tion of the components 'of force. The variations of the 

declination d& and dip di, are linear functions of dX 3 dY 

(i) (i) 

and dZ 3 and therefore of dA and dB , by virtue of 

n n J 

the known formulas:* 


Bd& =• — * sin &.dX + cos &.dY 3 

Bdi = — cos i sin i cos &.dX — cos i sin i sin S.dY. 

+ cos '^ m i.dZ 3 

We obtain these formulas and the following differential 
formulas, by differentiating the system: X = S cos i cos 6, 

Y = S cos i sin 6, Z = S sin i 3 or the system which we 
deduced by solvina 

6 and i (J? - S cos i ) with respect to B. 
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where R stands for the horizontal component of the 
terrestrial magnetic force. These are the formulas 
used by Erman and Petersen to reduce the observations 
to a common date in order to revise the Gaussian cal- 
culations they performed and which we will return to 
later in this report. 


Unfortunately, this method becomes unfeasible when 
a large number of terms are taken into consideration, 
because of the considerable amount of work required to 
solve the conditions equations. To clarify this, 

/ A ) 

let us assume that we have to determine coefficients 

(i) th 

arid B up to and including the terms of the i order* 

we would then have to solve a system of equations with 

i(i+2) unknowns. Moreover, the resolution of a large 

number of equations with 25, 35, 48 ... unknowns, using the 

least squares method, is certainly the typ j of task to 

discourage the most enterprising calculator, especially 

if he must repeat this process 20 times 


(i) (i) 

the variations of A\ and Zr 

n n 


ducts S'jj .3 


We might consider expanding d 5 and di directly as a 
Ci) (i) 

function of dA and dB , and I admit that this is ac- 
nn 

tually the method that I tried first in order to obtain 

We show that the pro- 

dS a 2 2 .dS 0 Z .di . , , _ 

S oos <5 qos may he expanded from 

the spherical functions. Indeed, these products are made 

j y 

up of a series of terms having the expression X^ J or 

X .X . ^ X representing any one of the quanti- ^ 

* c dt 

ties Xj Jj Z: these will be the expansion products accord- 
(i ) 

ing to the J' ' functions, and these products themselves 


will be expanded into a series of J 


(i) 


n 


The coefficients 


of these expansions will be made up of terms from one of 

- • .u> .w*®* 


the expressions A. 


A. 


l P 


The coeffi- 


•n dt J n m dt 

cients of these expansions may be expressed by set inte- 
grals, and since we therefore have enough data to solve 


the problem, their values will be calculated independently 
of each other by simple quadratures. Moreover, as we are 
assuming that these polynomials are known, in order to de- 
the v al ue S of ana „<<> . we wi U have to soXve 

a system of 2i(2i+2) or 3i(3i+2) equations for i(i+2) un- 
knowns, and we will see the same difficulty appear again in 
a new form. 


It was therefore necessary to use a simpler and more 
direct method. The method which we use consistently is 
a method of approximation, which consists of calculating 

accord- 


/ * j / * i 

the variations of the arbitraries dA\ %J and dB^i 

n, n 


ing to ordinary formulas which give the X , Y and 2 com- 

(i) (i) 

ponents as a function of A and B . This is achieved 

n n 

by adopting the values resulting from the first terms of 
th* 3 , formulas below as the approximative values of dY 
dZ: 


dY - 2? cos <5. <56 + sin6 dH 3 

dZ coA’ di + tan 9* 

the first terms are much larger than the second ones, if 
6 and i are small. We therefore disregard the variations 
of intensity which were not revealed by earlier observa- 
tions. 

A second approximation, in which. 2? is variable Cwe 

may use the value resulting from the first calculation) , 

(i) (i) 

will then give a more approximate value of A n and B n , 

and by continuing in this way we may extend the approxi- 
mation as far as we want to. 

This method does not perhaps always offer the most 
reliable means of obtaining a calculation of the varia- 
tions of the arbitraries, particularly with respect to 


the dip. However, from the standpoint of application, 
it offers great simplicity. 


2 . 


In the first approximation, we limited ourselves to 
the terms of the first two orders and therefore set: 

by noting that for a magnetic body « o. 

The complete formulas used to determine the unknowns 
will be given later in the text (no. 6). Let us simply 
state here that we have followed the Gaussian method to 
perform our calculations, by evaluating the sin >inr and 
cos met coefficients in the expressions representing the 
forces by means of mechanical quadratures, after which 
the unknowns were divided into several groups in order to 
make it much more convenient to determine them. 

The calculations were extended successively to earl- 
ier and earlier periods in history, starting with the year 
1829. : for. the declination, from 1829 to 1787, from 1787 to 
1700, then from 1700 to 1600; for the dip, from 1829 to 
1787 and from 1787 to 1700. 

Having obtained, with the accessible data, the values 
of the declination and the dip at the intersecting points 
of nine equidistant parallels, between the northern 
and southern latitude of 60°*, with 12 meridians in the 
longitudes 0°, 30°, 60° East of Greenwich, the differences 
i? oosS ,.dS or were calculated. Having no knowl 

edge at all of the intensity variations, we were forced to 


*For 16 0.0, and for the dip in ,1787, observations of 
the parallel 60°south were lacking; for the dip of 1700 
the parallels of 60° north and south were both lacking. 
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in this first approximation, the 18 29 values for if, a. 

which inevitably would introduce errors that the sub- 
sequent approximations had to correct. These differ- 
ences were then assimilated with dl and dZ, and sub- 
stituted in the formulas of no. 6, reduced to their 
first terms. 


After calculating the values of coefficients M and 

N using the method just described, I then found the var- 

iations d'A^ 1 and dB„- for the selected time intervals, 
yi n 

which when successively added to the intial values for 


(i) 


(il 


1829, gave the values for A and B for the selected 

VJ VJ 

dates . 


n 


n 


One important condition required for the validity 

of our method of approximation is that we can disregard 

(i) (i) 

the squares and the products of variations dA n and B^ 
This condition is far from being fulfilled for these 
long intervals of hundreds of years. We see to what ex- 
tent the causes combine to distort the results, and we 
should expect no more than a very rough approximation of 
this first test. 


3. 


Our numerical calculations are based mainly on the 
magnetic charts accompanying Hansteen’s work, Nn ter such - 
ungen uber den Magnetismus der Erde, published in Christ- 
iania in 1819. The very incomplete charts of this author 
for 16 QO. and 1700 were completed/ in the Pacific, with 
information at our disposal. 

For the 1600 chart, we took, advantage of Tasman's 
observations, made in 1642 and 1643, near New Holland. 
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This included W. van Schouten's observation southeast of 
the Marquesas/ reported by Hansteon. Another was Candish’s 
observation/ cited by Kircher in the appendix to the 2nd 
edition of his work Magnes / sive de Arte magnetica (Col . 

Agripp. 1643) . Finally/ we used another observation in the 3rd 
table of oceanic observations. Along the northern coast 
of America, my attention was drawn to a compass-card on a 
Wytfliet’s chart, in Fascimileatlas by Nordenskiold and 
which I identified after noticing an error in a sign. 

For 1700, we have observations by Wooden Rogers dur- 
ing a trip from California to the Philippine Islands in 
1710, and a few observations east of New Holland. 

The chart of Wilcke'* equal inclination lines is 
found in volume XXIX in M£molres de l'Acaddmie de Stock - 
holm, 1768 . The isoclinal lines of the North Atlantic 
are plotted here in full outlines according to Ekeberg 
and Lacaille’s observations and in broken lines accord- 
ing to P. Fevillde's observations in 1710-1711. I followed 
these-. In the Indian Ocean, Wilcke had plotted the iso- 
clinal lines according to Ekeberg ' s observations in 1766. 
However/ in these regions, I preferred earlier observa- 
tions made by Cunningham in 1700. Finally, in the Pa- 
cific Ocean, a very precarious interpolation gave me very 
inadequate values of the inclination. I have to admit 
that this part of the curves is based more on conjectures 
than on reliable facts. 


4. 


I believe I have given a full account of the numbers 
found, as provided directly by the first calculation/ al- 
though they have no other value than to show an example of 
how one can pass from total ignorance to a more and more 


rigorous knowledge of the truth; 


on/or.u p/igb j s 

OF POOR QUALITY 


The numerical values of the 8 constants of the first 
two orders, obtained by the first calculation, are shown 
below (the second lines corresponding to 1700 and 1787 
refer to inclination observations) : 



Years 

< 

b\' 


i6co 


—0,018058 4* 0.023456 


1700 

— 

— 0,032392 1-0.031618 

* 

» 

+ 0,288411 

—0,049858 4-0.080636 


1784 

4- 0,341369 

—0051882 4-0.032882 


1787 

— 

—0049249 * 0028367 


> 

+ O.325972 

—0.042256 4-0.051105 


1829 

+ 0.320074 

-0,060109 4- 0,0^353 


1830 

+ 0-323477 

—0062456 4-0.031106 


1858 

4- 0,324627 

—0,057268 4- 0.023641 


1880 

+ O.333923 

— 0,061920 4-0.027636 


t88; 

4- 0.315720 

— 0,060258 4. 0.024814 

Years 

Bo' 


U\ A n n i 

)6oo 

— 

1 

0 

0 

M 

,co 

U 

— 0042494 4 - 0.026955 4- e.016715 

1700 

— ■ 

—0.012273 

- 0.O36263 4- 0.005250 4- 0.014823 

) 

+ 0,023114 

4- 0,006626 

-0006325 —0.000898 4- 0.022994 

1784 

+ 0.004694 

— 0,018639 

-0030099 — 0.013829 4- 0.0:6553 

1787 

— 

— 0,015681 

— 0037125 —0,0056844-0.007518 

9 

— 0.014242 

— O.OU230 

— 0,03 1 8 5 5 — 0.0 1 1 063 4* 0.008777 

1829 

4- 0.001 210 

4- 0,000720 

— c 044537 — 0.012637 4 0,001249 

1830 

— 0.007708 

— 0,002107 

— 5050635 — 0,013631 4-0x0017? 

1858 

— 0.007190 

— 0.002226 

— 5.050827 —0.015093 — 0x01337 

1880 

— 0,005507 

— 0,003805 

— 5051930 — 0.012433 — 0.004284 

1885 

4 - 0,007906 

+ 0,012999 

— 5 C49798 — O.OI 2604 — 0,005667 


We have inserted in this table the results of earlier 
calculations, performed before ours, and in particular by 
Gauss for 1830 , by Erman and Petersen for 1829 , by Mr. von 
Quintus Xcilius, for 1880 , and finally by Neumayer and 
Petersen, for 1885 . We have included Erman' s and Peter- 
sen's values for 1784 and 1858 , deduced from the parabolic 


formulas they gave in their report on page 11. 


Here are the numbers which will serve as the basis 

for our subsequent calculations. 

For those who are interested in referring to the 

original sources/ the relevant bibliographical informa- 
tion is provided below: 

-Gauss, Allgemelne Theorie des Erdmagnetismus (Gauss 
Werke, t. V; Resultate aus den Beobacthungen des mag- 
netischen Vereins; 1838). 

-Gauss & Weber, Atlas des Erdmagnetismus nach den Ele - 
menten der Theorie entworfen , Leipzig 1840. 

-Goldschmidt, Vergleichung magnetischer Beobachtungen /9 

nach den Blementen der Theorie , (Resultate aus den 
Beobachtungen des magnetischen Vereins; 1840) . 

-H. Petersen, Vergleichung der gaussischen Theorie des 

ne ti smu s rait A. Ermans magnetischen Beobachtungen, 
fk^tronoraische Nachrichten, no.s 450, 452, 454). 

-Erman & Petersen, Die Grundlagen der gaussischen Theorie 
und die Erscheinungen des Erdmagnetismus im Jahre 1829 , 

Berlin 1874. 

-G.Von Quintus Icilius, Der magnetische Zustand der Erde , 

(Aus dem Archiv des Deutschen Seewarte, 4th year, 1881, 
n . 2 ) . 

-Neumayer, Atlas des Erdmagnetismus , in the 4th section 
of the new edition of the Berghaus Atlas on physics, 

Gotha 1891. 

-Ad. Schmidt, Mitteilungen liber eine neue Berechnung 
des i-irdmagnetischen Potentials . (Abhandlungen der konigl . 
bayrischen Akademie der Wissenschaf ten, Vol. XIX, 1st part, 

Munich 18 95. ) (For translation see N ASA- TM- 77329) . 
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These values contain all sorts of errors due to 


the imprecision of the charts and to the influence of 
the terms omitted in the expansion of the A and b]^ 
values according to the powers of time. In order to 
reduce the detrimental influence of these causes of er- 
rors as far as possible, it is necessary to associate 
the numbers found in a formula expressing the law of 


variation for constants A 


(i). 
n 


and B 


<i) 


n 


We know nothing <2 about the form of quanti- 

(i) (i) ' 

ties A and B taken into consideration as a function 
n n 

of time. It seems natural at first to expand them into 
a series according to the powers of time, and to stop 
the expansion at a power value that is not too high. 
However, a moment's reflection will tell us that these 
quantities cannot contain secular terms per se in which 
time is a factor. Actually, this would imply that the 
quantity y^’ 2 + and consequently the corresponding 

magnetic moment could become as large as we want it 
to become, which is impossible. 


It is also not very probable that they will get 
close to a limit, as we may assume that the earth's mag- 
netisation dates back far enough for a steady state to 
be reached in the variations. 


The only probable hypothesis is that they do not 

really contain any periodic terras. One of the simplest 

hypotheses that we can make in this regard is that each 

(i) (i ) 

pair of transcendental quantities A n and % n has cycle 
which last the same period of time. It remains to be 
seen that this hypothesis is verified by experience. 


/•-,*) (i ) 

Let us replace the quantities A‘ n ' and F._ i(i+2) 


n 


12 



ORIGINAL* PAGE r$ 
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by new transcendental quantities a nd ^ related 

to the first ones by equations having the expression: 


tf'K 

n 


41 r 


7 =* , 


E 


rp ,1 !■' _ ) 


where E represents the base of the Napierian logarithms. 

will be the module or the magnitude of the magnetic 
moment )/ A™* + 2?'* , angle will distinguish its dir- 
ection in space. n 

(i ) , 

Once this is established, the part of I depend- /10 

end upon angle «nr, will be expressed: 

0 — • “? COS II (O + . 


We will look for the law of variation for the new 
(i) (i) 

constants a and ft • 
n 

( 'V ) ( ) 

The numerical values of the constants a and 

derived from the values of A [% and B' found in the 

n n 

preceding number are shown below (the g may always be 

( 'I* ) 

selected in such a manner that the values of a will be 

n 

essentially positive, if n > o) : 


U) 


and 




The values of the coefficients a 

( 1 ) ( 2 ) 0 

identical to those of B and B given earlier. 

o o ^ 


are 


These numbers show more clearly the law of variation 

(i) (i) 

of the arbitraries and 6^ and we already guess that 

this law is very simple. Actually, the variations seem 

(2 ) 

to be proportional to time. It is the function Y 0 in /II 

" “ ^ (£} 

particular which gives us some evidence: coefficient ^ 

*We may assume that this angle is equal to 133.67, if /10 

we give the corresponding coefficient a-£ 2 ^ a negative sign. 

The value of given in the next column (p.14) refers 

to the value of g.J^i minus 180°. 
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remains fairly constant, while angle has shifted /ll. 

nearly a quarter of a circle over the past three centur- 
ies and this me -ms that if it keeps up at this rate it 
will take a total of 1200 years to move around in a full 
circle . 


Years 

_(i> 

“o 

W P '- P > 

ifoo 

— 

+ 0029602 + 37°.6o — 7’. 1 8 

1700 

— 

4 0.045970 4 44 .80— i .34 

> 

4 0.2884 II 

4 0,094805 + 31 .73 4 11 .73 

1784 

+ 0,341369 

-r O.C61424 + 57 .63 . — 3 .20 

1787 

— 

4 0056836 + 60 .05 — 5 .25j 

> 

+ a 3 2 S 97 2 

4 0.0663 13 + 39 -58 + 15 -24 

1829 

4 0.320074 

- 0.066461 4 64 .75 — 4 .45 

1830 

+ °- 3*3477 

4 0069773 4 63 .52 — 3 .09 

1858 

+ 0.324627 

4 0061956 4 67 .57 — 3 48 

1880 

+ <* 3339 2 3 

4 - 0067808 4 65 .95 4 I .01 

OO 

OO 

4 - 0.315720 

ra \ 

4 0065169 4 67.62 — O.OI 


Years 

Ji ) 

K l) 

<a> vs> 

, «i P\ 

pc — P 0 

PT Pc 

@0 

1600 

— 

4 0.051084 4 I46'.2S 

— 1 7 ~‘ 7 ° 

4- 0.031724 — zg\io — 

i ’,76 

1700 

— 

4 0.038284 4 161 .30 

-13.15 

4- 0.014789 — 10 .40 4 

10. 11 

7 

4 0,02 3114 

4 0.009160 4 313 .67 *■ 

4- 14 - 4 ® 

4-0.0230114- 1.12 — 

1. 41 

1784 

4 0.004694 

4 0.035404 4 148 -23 

4- 16.37 

4- 0.015303 + 32.33 — 

6.93 

K 

00 

K 

— 

4 0.040301 4 157 .10 

4-8.08 

4- 0.C09425 + 18 .55 4 

7.76 

> 

— 0.014242 

- 1 - °-°33777 + .5S 

4-4.60 

4- 0.014122 4 Z5 ,79 4 

0.52 

1829 

4 0.001210 

4 0,044543 + 1S0 ,94 

- 7-53 

4 O.OI2699 4 42.18 — 

3.02 

O 

OO 

— O.OO7708 

4 0.050677 4 177 .62 

— 4. 02 

4 0.013632 4 44,64 — 

S .i8 

1858 

— 0^190 

4 0.050876 4 - 177 .50 

4 I .58 

4 0.015152 4 47-54 + 

0.48 

1880 

— 0.005507 

4 - 0.052068 4 175 .82 

~ 7 -S 6 

40.013151 4 54.51 + 

0 .25 

1885 

4 0.007906 

4- 0,051468 4- 19465 

— IO .26 

4 0.013819 4 57.11 — 

0.83 


This leads us to represent the variations of the 
constants using linear formulas of time. We have found 
that if we give each determination the same value: 


I 
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«0 S3 + 0.322J92 + i‘0.10^, 

a \ l> S 3 + 0.064528 + 49,46 1 , 
rtf” s= —0.000389“ 68.21 t, 
«**’ a 4. O.O453II + U3.36 t , 
«»” » +O.OI3577- 57 - 3 <>*f 


/ 9 ( , u » 5 o’, 30 + o". 1305 1 , 

a 173.40 + 0.1957^ 

= 39.15 + o* 3 °S^ ** 


Time is counted as of 1829, and the variable terms of 
the modules are expressed in units of the sixth decimal 
order. 

( A ) 

Coefficients a„ vary very slowly with time, and 
it is permissible to assume, within the limits of obser- 
vation errors, that the time factor is strictly zero. 

The very large discrepancy between the numbers found do 

(i) 

not enable us to determine whether the a modules 

n 

still contain periodical terms. However, for theoretical 

reasons which we will give later (no* 15), it is probable 

(i) 

that the periodical inequalities which may exist in a 

Yl 

are of the same order of magnitude as the periodical in- 

( *£• ) 

equalities of angles $ n As the observations have 

still not revealed any such thing, we have decided to set: 

(i) 

a = const. , 

n ' 


and therefore: 


Y, 


Ci). 


n 


and m 


a 2 


n 


. a CU _ JU , m U) + 

K ~ Y n + n n 
being two new constants. 


/ • \ 

By determining the or 's according to the six most 

n 

accurate modern determinations, and by excluding pre- 
vious observations which mainly produce the term in t, ' 
we have finally set:* 


*A comparison of these values with the small table below 
shows that the declination chart for 1600, and that of 
the inclination for 1700 are highly inaccurate, as we 
could naturally assume. 
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«0 1 » + 0-3^533 . 
n\ U a •{* 0.065432 , 
«r « — O.OOF099, 
rt < * > = + O.O47J06 , 
< tf } a + 0,013959! 


fl [ x) ** fo’,30 + o’. 1 305 t , 

PT ™ >73-40 + 0.1957#, 
P \" = 39.15 + 0,3058 


These will serve as the base values for our second 
approximation, in the calculation of the A ^ 'e & 's 

as of the year 1538 to and including 1800. 


6 . 


Once we have found our first approximation, it will 
be extended as far as possible. The same formulas are 
applicable to the case already processed, by letting the 
characteristic d no longer be the difference with respect 
to time, but the correction which must be made of the ap- 
proximations <cbf the uiiknown values. 

To these rectifying calculations, we have extended the 
series of the potential up to and including the terras of 
the fourth order, which are indispensable for a correct 
representation of the facts. 

These complete formulas which were used to calculate 

the DA^^'s and the dB^^'s according to the d.Xs 3 dYs and 
n n 3 

dZs will serve at this point. 

If we set: 


Py = (i-u-)2 


J Mi 


dfi 1 * 


the general term of the function for the magnetic forces 
will be (no. 1) : 
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M' 

yi+i 


J?/. (J./ sin n& -f P/ cos ntr) , 


In the expressions of the components, the terms 
corresponding to any given point on the terrestrial 
surface will be: 


1 * ^ ,,ti3r + COS >m ) f 

~ cc ' 5 nvx ~ ■®i 0 sin n ™) ? 

(* + 0 <*?• (-4° sin m + 1 S|'» cos nm) . 


If we call ■ M^ } the coefficients of 

cos nw and sin nw in the expansion of X s Y and Z accord- 
ing to the multiples of nr/ the equations which determine L 


M n * etc., 

will be: 




i_ 

II 


x.= 

V n pd> 

Jm \l-u *- “ ' 

-e 



K= 

V }i 73(0 



By focussing on the terms up to and including the 
fourth order, we will find in particular the following 
formulas: 
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\T^i ? . + 2,« VT^-T 5 . J3? + (3ft 2 - “ « 2 • B f + (4f* 3 - VP) V' « - f‘ 2 • ^ » 

2 fi . Bf + 3 (f< 2 ~ J) - *? + 4 (f* 3 “ Isf 1 ) • + 5 (A 4 ~ tM 2 + A) • -Bo* > 

— li . 4> — (2ft 5 — i) . -4? — (3ft 3 — Vf‘) • A T (4P* — VP * + f) • ^ 4> * 

„4” -ft.4 2> -(M 2 -jt)-4 8) -^-ffi).4 4; , 

2 ^=^.4 l) + 3f<y7=^ . 4 2) + 4(f< 2 ~0 + 5(A 3 -?P)V / i-P a -4 4) i 

= — ft, JBf — (2ft 5 — i) .3® - (3ft 3 - ¥f*) • ■Bf — (4ft 4 — Vp 2 + $) • Bf > 

= 4‘> + fi. jef + (m 2 - i) • Bf + (p 3 — ?p) • Bf > 

= oyTZ^l . Bf + . Bf + 4 O' 2 — i) VT^ 2 . Bf + 5(ft 3 - ?f0 V'^-ft 2 • Bf ? 

= - 2 f(»'T=^ . 4 2> - (3M 2 - i)v"T=7 2 • 4 3) - (4^ 3 ” M ^ t=: p 2 • A f * 

=s _ 2 yl[tr^.4 2 > -wYT^F.A? 

_ 3 (I — ft 5 ) . 4 2> -f 4ft (i — ft 5 ) . /if + 5 (ft 2 — 4) 0 — P 2 ) • 4 4) ’ 

= _ 2 ftvr— 7T 2 . i*f - on* 2 - 0 • B f - (4,« 3 - Yp) v'r-^ • Bf • 

_ 2)1—0* . Bf + 3ft V^-T 3 ■ Bf + 2(fi 2 — I) — P* ■ Bf » 

== ‘ 3(1 -ft 5 ) . JBf + 4.“ (I -P 2 ) • Bf + 5(ft 2 - 0 (I -M 2 ) • Bf i 

« — 3 ft (i — ft 5 ) . 4 3) - (4ft 3 - i) (i - ft 2 ) - 4 4> * 

= _ 3 ( I — ft 5 ) . 4 3) — 3ft (r - ft 2 ) • 4° > 

= 4 (VT4T 5 ) 3 . 4 S) + Sft (/F=^) 3 . 4 4> ; 

= - 3ft (I - ft 2 ) • 4” - (4ft 2 - I) (I “ ft 2 ) . -Bf , 

= 3(1 — ft 5 ) . Bf + 3ft (I - ft 2 ) . Bf , 

= 4 ( yr=^) 3 . tff + sft^i^y.-Bf? 

= - 4 ;<yr=^) 8 .4 4) , 

- -4(vr=^) 3 -4 > » 

- 5 (i -ft 2 ) 2 . 4; 

= -4fta / ^) 3 -4 > » 

= 5 <i-u 5 )Mjf. 
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My first test made me aware that the charts I used 
were much too inaccurate, the large gaps left behind 
by ancient observations were filled more by bold conjec- 
tures than by verified facts. Specifically, when I mis- 
interpreted a compass-card on Wytfliet's chart, I gave 
the Pacific coast of northern America a western declina- 
tion rather than an eastern declination and thus made a 
serious error which calculations quickly corrected. 

Since then, X found several ancient observations 
which enabled me to construct new charts for the XVlth 
and XVI I th centruies, and specifically for the years 1538, 
1572, 1600, 1642, 1676 and 1700. The charts for 1600 and 
1700 encompassed a complete region of the earth, including 
all longitudes. For the other periods, the dates needed 
in the Pacific were lacking: they encompassed only the re= 
gion of the earth between longitude 90° West and longitude 
150° East. 

We will not reproduce these charts containing too many 
errors, we will simply show later in the text the results 
of a comparison of the chart for 1600 with theory. 

We will succinctly indicate the main observations 
which served to construct the new magnetic charts. 

The first of our charts is based mainly on the obser- 
vations made by D. Joao de Castro during his trip from Lis- 
bon to Goa in 1538 and to the Red Sea in 1541 ( Roteiro de 
Lisboa a Goa , published by Joao de Andrade Corvo in Lisbon 
in 1882; Prlmeiro Roteiro de Costa da India , ed. Diogo 
Kopke, Porto 1843; Roteiro ao mar Roxo no anno de 1541 , 
Paris 1833) . The second group of charts are based on the 
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observations of the famous Christopher Columbus during 
his third voyage. The charts are also based on various 
notes from Sebastian Cabot's voyages and on Magellan's 
first voyage around the world (Pigafetta, Navigazione in 
torno al globo , reprinted in Raccolta di domuneti etc . 
dalla R. commissione colombiana , Part V, Vol. Ill, or in 
Atti della Socita ligure di Storia patrie , Genova, t. XV, 
p. 31). Finally, our charts are based on a certain number 
of data on the cities of Europe given by Tanstetter, 

Pierre Apien, George Hartmann, Joachim Rheticus, Mercator, 
Bellarmatus, Orance Fin6, Cosimo BartoXi, etc., and others 
who we could mention, but which would make our list too 
long if we gave them all here. 

The major gaps left by these observations in vast 
regions can be filled only by a very precarious process. 
For this purpose, I used ancient mariner's charts. We 
know that the ancient navigators were accustomed to con- 
structing their mariners charts according to compass indi- 
cations without being concerned about variations in the 
magnetic needle. By comparing the ancient mariners charts 
with modern charts, we can deduce the declination of the 
magnet. Among these charts, we will mention those which 
existed in Harrisse's work ( Discovery of the North Ameri - 
ca , Paris 1892): the maps of la Cosa, Cantino, Vesconte 
Maggiolo, Estevam Gomez, Canerio, and a planisphere kept 
in the library of the king of Italy in Turin, dating back 
to the very f irst years of the XVIth century. We may also 
consult fruitfully the Facsimileatlas by Nordenskiold, 
passim . and especially the charts of Bernardus Sylvanus, 
Martin Behaim's globe, etc. 

The following charts are based on somewhat more re- 
liable observations. For 1572, we have the observations 
of Vincente Rodrigues, published by Andrade Corvo (loo. 
ait,), according to a manuscript at the royal library of 
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Madrid. We also borrowed from a list of declinations 
extracted from Art de trouver les ports iThe Art Of Find- 
ing Ports] by Simon Stevinus. II had under ray eyes Edw. 
Wright's translation of the third edition of his Treatise 
Certain Errors in Navigation Detected and Corrected (Lon- 
don 1657)]. This Table was often reproduced by the major 
compilers of magnetic observations of the XVII th century, 
Gilbert, Kircher and Riccioli. We will also join a 
few facts on Adrian Metius given to us by P. Kircher in 
his famous work ( Magnes, sive de Arte magnetlca ; p. 429 
edition of Rome 1641).. A certain number of European ob- 
servations taken here and there will complete the chart 
for this date. 

In regard to the next chart, fori 1600, we tried to 
make it as complete as possible, and were essentially 
successful in our attempt, thanks to the observations in 
the South Sea available in Dudley's sea Atlas, called 
Arcano del Mare , Florence 1646-16.47, and particularly to 
those of the southern hemisphere made by van Schouten and 
Lemaire in 1616. In the Old World, we have borrowed from 
observations made by Gaspar Remao in 1598 and published 
by Andrade Corvo as an appendix to D. Joao de Castro's 
voyage. We also borrowed observations by Thomas Candish, 
Hall, Davis and Linschot, in P. Kircher ' s Magnes (appendix 
to the second edition. Col. Agripp. 1643)., an a large num- 
ber of observations by Middleton, Daunton, Hippon, Saris, 
Castleton, Marlowe, etc., in the 3rd volume by Purchas 
I His Pilgrims (".London 1625);, quoted from Hansteen, Unter - 
suchungen uber den Magnetismus der Brde , Christiana, 1819.]. 

For the other two charts of the XVIIth century, we 
will find various observations reported in general works 
of the period, by Kircher, Riccioli, Fournier, Semmlerus. 

For the first half of the XVIIth century, we will 
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cite observations made by Fournier and Guerard tG. Four- 
nier, Hydrographie {Hydrography;] , Paris 16431. Guerar- 
dus Diepensis ' observations are also cited in Riccioli 
( Geograph ia et Hydrographia reformats , yenetiis 1672)., 
who also copied all of his predecessors, Kircher, Dud- 
laeus, Stenvinus, Janssonius, etc. We then offer two 
important observations by Olearius in the Caspian Sea 
(.Dells le, His to ire de l'Academie , I history of the AcademyJ 
1721.) Tasman’s observations date back to the same per- 
iod. For the North. Atlantic, consult Bichard Norwood, 

The Sea -mans Practice , London 16'59.. In Europe, we bor- 
rowed very numerous observations gathered by the diligent 
P. Athanase Kircher, and published by him in his work al- 
ready cited so often, as well as observations in the Med- 
iterranean basin reported by the same author. 

The 1676. chart is partially extracted from the same 
sources as the previous one, and especially from Christo— 
phorus Seramierus' book, Methodus inveniendae longitudinis 
. maritimae , Halae Magdeburgicae 1723. Let us also mention, 
for the end of the XVlTth century, Halley's Table in the 
Philosophical Transactions of London (volume XXII, 1683).; 
Leydeker' s observations in 1675 (Hansteen, Untersuchungen 
etc.).; and a host of French Jesuit observations made in 
France, in China and during the trip to China ( Anciennes 
Mdmoires de l'Academie , t. VI I, passim.)’ I Ancient Beports 
of the AcademyJ 

Our chart for the year 17 00. is, in all of its main 
aspects, a well-known reproduction of the Table by Moun- 
taine and Dodson in the Transactions philosophlqUes de 
1757 iPhilosophical Transactions of 1757J , p. 329. In 
the Pacific Ocean, the curves were plotted according to 
the voyage of Captain Wood Rogers 1710 (.T ransactions phil - 
osophiques, 17 21, p. 173)^, and in the northern part fol- 
lowing the '.voyage. of the ship Navire St Antoine in 1707, 
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based on an interesting description of this voyage 
preserved in a manuscript at the Royal Library of 
Stockholm and which I became aware of through Mr. 

Dahlgren. The observations of Bering at Kamchat- 
ka made during his voyage in 1728, are extracted from 
handwritten charts of Bering, also preserved in 
the Royal Library. Inside Siberia, we have only Gmelin's 
observations in 1735, and a few notes by Euler, of a 
fairly uncertain date, plus van Verden's observations in 
the Caspian sea (Histoire de l'Acaddmie 1718). For Eur- 
ope, consult Musschenbroek ( Dissertatio de Magnete , Viennae 
1754), which also reports a host of nautical observations 
for the turn of the XVIIIth century. By including a 
certain number of observations in the Hudson Bay in Mus- 
schenbroek and a few other American observations in Ch.-A 
Schott ( Secular variation of the magnetic declination in 
the United States ; Report of the tUS. Coast and Geodetic 
Survey , 1888, App. n. 7, 7th edition, Washington 1890), 
we believe we have indicated the most important sources, of- 
fering our chart enough background to give it a broader 
extension, something which was not possible in Halley's 
time. 


Here are the main sources from which we have extracted 
our information. We could not invite the reader to follow 
us in all the often tedious research details which we went 
through to establish the location or date of the observa- 
tions which ancient authors often described with vague de- 
tails. 

In conclusion, I' would like to point out that in the 
charts which I had decided to sketch, I did not use any 
of the numerous and precious observations recently dis- 
covered by Mr. van Beramelen, in the Dutch Record Office*, 
in fact I was unaware of his work when I set out to 

*See his inaugural thesis presented before the 
Science Faculty at Leyde, 1893. 


to construct my charts. Fortunately, there is no roason 
to be overly regretful about this unintentional omission, 
as X am convinced that the probable differences between 
his charts and mine are generally between 1° and 2°. This 
is the error range to be expected in ancient observations, 
and the errors which may be found in our conclusions are 
of little consequence. 

Besides the six new charts which we had undertaken to 
sketch, we referred to Hansteen's well-known Atlas where 
equal declination lines will be found for the years 1710, 
1720, 1730, 1744 , 1756, 1770, 1780, 1787 and 1800, together 
with Wilcke's inclination chart for 1700-1760. 

Finally, three new determinations of the potential 
for the XIXth century will complete the series of unknown 
values which will serve as a base for our studies. These 
determinations for the years 1820, 1840, and 1860, are 
based on the following data: 

1. Hansteen's declination chart ( K. Danske Videnskabs 
Selskabs skrifter, 5th series, t. IV; 1855) , and the incli- 
nation chart by the same author for 1827 ( Annales de Pog- 
gendorf f , t. XCVII; 1830), attached to Duperrey's total 
intensity chart, based on all intensity observations made 
from 1784 to 1828 ( Bergahus' physikalischer Atla s, 2nd 
edition. ) 


2. Sabine's total intensity charts for 1837-1840, 
attached to the declination and inclination charts by the 
same author for 1840. ’The chart of isodynamic lines ap- 
peared first in the Report of the British Association , 
1837. The declination chart is found in Johnston's Phys - 
ical Atlas . All of these charts are reproduced in F. 
Walker, Terrestrial and cosmlcal magnetism, Cambridge 
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3. Finally, Evan's charts of the three elements of 
terrestrial magnetism for 1860 (Smith, Admiralty manual 
of the Deviations of the *, 3rd edition, London 1869). 

8 . 


/ • \ f • | 

The calculations of the Ai 1 8 and B„'b based on 

n n 

these data were performed according to the same procedure 
used for the first approximation, by utilizing the complete 
formulas of no. 6, and by following the same calculating 
rules. These will not be repeated here. 

However, the incomplete charts, where the Pacific was 
missing, required another method for calculating the coef- 
ficients Lj M and N which could no longer be obtained using 
ordinary quadratures. We have obtained the values of the 
nine unknowns by solving the system of nine equations pro- 
vided by the points of intersection of the nine meridians 
with each parallel. 

These equations have consta t coefficients, and it 
is possible to benefit from this situation to abbreviate 
the calculations. Let us call the numerical 

values of the second members of the condition equations 
for o” = 0°, 30°, ...150°, 270°, 300°, 330°, used in the 
order just enumerated. We will immediately express the 
coefficents of these equations. By solving them, we ob- 
tain the following result where the coefficients are log- 
arithms: 


♦illegible . 


mmmi pass eg 
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L 9 0.3336?^ + 0,3960 n a — 0.3336 »j + 0.271 m«— 0.0324 n» + 1.794* ««+ > -7939 «7 — 0.0326 n* + 0.27 1 1 w a 


//, 0 + 0,5721 

-O.572I 

+ O.572I 

—04690 

+ 0.2710 

— 0.0327 

-78971 

+ 0,2711 

—04242 

L \ S* + 0,3336 

-O.3335 

+ 0.3335 

- 0.1357 

4* 0.O324 

—74612 

—78970 

+ 0,0325 

—0.271 1 

l t a —0.1356 

+ 0,27 JO 

— O.27IO 

+ 0,1356 

— 0,0730 

+ 78346 

— <$b 

—78346 

+ 0.0732 

L j 0 —0,4690 

+ O.5O96 

— 0,4243 

+ 0.3336 

— 0,1679 

+ 75966 

+ 7.8970 

—0,2458 

4 - 0.3544 

ijK — OO 

— OO 

— 00 

-74603 

+ 76990 

—74603 

+ 74606 

-76991 

+ 1,4606 

+ 0,3335 

- 0.3335 

+ 0.3- *- 

—0.2711 

+ 00325 

—7,4607 

—7.4618 

+ 0,0326 

—0,2710 

Li => + T.897O 

-7.7939 

+ 7,46*. 

— OO 

-“0237 

+ "2,7846 

— 7 o 86 o 

+ 75962 

—7.8346 

l .\ a — 7,8970 

4 - O.0324 

— 0,0324 

+ T 8970 

-7,5959 

+ 70241 

— OO 

— 70249 

+ 7,5960 


Similar equations determine the values of M n > M^j 


Besides the four new charts, this method was used 
again for Hansteen's charts 1720, 1730 and 1756 and for 
Wilcke's chart for 1700-1760. The rest of the calcula- 
tions were performed as in the preceding case. 

The results of these calculations, in the form direc- 
tly obtained, are included in the following tables (p. 18 
19).. These are fundamental data, directly deduced from 
the observations, which will serve as the base for all of 
our subsequent reasonings. 


We will make the following remark in regard to these 
figures. Wilcke's inclination chart for 1700-1760 can 
hardly be used to determine the terms of the first two 
orders. This is because the insertion of higher order 
terras distorts the result. The numbers given below are 
those resulting from this assumption and where 3rd and 
4th order terms are omitted.* 


*We will cite, for the sake of curiosity, the numbers 
which we had obtained first, by using the complete for- 
mulas, for the independent functions of longitude: 


Dates. rf " J 3 ( 0 4) 

1700 4 - 0.1*978* + 0.143782 

*760 + 0.38706a + 0.244060 


1/ 

O 0 "0 

— 0.321506 +0.153582 

+ 0,055931 + 0.6238*4 


These results , which are far from being accurate, show that 
we should. be very cautious about giving expansions a broad- 
er extension than the observations can substantiate. 


a ) 

Next, we calculated the module a' and the argument 

B u ' ; of each term, using the formulas of no. 5. The re- 
ft 

suits of this work, are shown on the next four pages. 
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In order to examine the hypothesis offered above on 

the law of variation of the transcendental quantities A n 

and we have calculated, as before, the coefficient 

(i) n 


a ' " ' and the argument 6 
ft ft 


(i) 


of each term. The results of 
this work are shown on the next four pages. 


These tables enable us to account for the variation 

of the unknowns with time. The linear variation of the 
/ »• \ 

arguments 8^ are obviously confirmed for the terras of 
the first two orders. In regard to higher order terms, 
the results are less obvious. However, as a general rule, 
we can affirm that the results do confirm the hypothesis 
used as the base for the second approximation, so that we 
may now set: 


= const.; 3^ 


ft 


m n ^ 


where t represents the number of years as of 1800. 

(i) 

If we refer to the meaning of the quantities and 

6^', the law of variation of the earth's magnetization 

ft 

can be expressed 'succinctly in the following words: 


"The modules of each term of the magnetic force func- 
tion remain constant, while their arguments increase pro- 
portionally with time." 

(/t) 

The numerical values of the initial coefficients ot M 

( * 1 ( ¥ b) n 

and arguments and their yearly motions m n are shown 
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years 

il'o 1 ’ 

,(1I H ‘I 

•*1 "1 

J ? 

AT 

>538 

— 

—.0,044672 + 0.0562-7 

— 

—0,068463 —0.005050 

1572 

— 

— O.06254I *f 0,056953 

— 

+ 0,000331 —0.045492 

j6oo 

— 

— C.O21374 +0,032864 

— 

—0.033729 —0,030416 

1642 

— 

— 0,093116 +O.06535O 

— 

—0.031050 —0.028134 

167 6 

~ 

— 0,075413 +O.068&6 

— 

—0024094 —0.026200 

1700 

— 

— OO2862O + O.O32-24 

— 

—0,022398 —0,033977 

id, 

+ 0.2882 J I 

— 0.059880 + 0,072848 

-0,063404 

— 0.026486 + 0.0833 1 3 

1710 

— 

— 0,031975 +O.032J25 

— 

—0.014468 —0,027477 

1720 

— 

—0,054521 +0,040967 

— 

— 0,037366 —0,028640 

*730 

— 

— 0045669 + O.033366 


— 0.042686 — 0.025Z16 

‘744 

— 

— 0,041497 +O.O32746 

— 

—0.025255 —0.040958 

1756 

_ • 

— 0.036876 + 0,0249 1 3 

— 

— 0.09418Z + 0,00525? 

1760 

+ 0.347964 

+ 0,004396 + 0.068x3 

+ OO24643 

—0,109768 —0.094125 

1770 

_ 

— 0,047369 +0031091 

— 

—0,014586 —0,037508 

1780 

+ 0,319216 

.— 0.044853 +0.047254 

— o,oc6ioo 

—0.008785 — 0.036419 

1787 

— 

— 0,050481 + 0,03101; 

— 

—0.015738 —0,038748 

iSco 

— 

— 0,052754 +0,034429 

— 

— o.oiozco — 0.04J947 

1820 

r 0.333810 

— 0,062973 +0.026552 

+ 0,005092 

—0.000091 —0,030416 

1829 

4 - 0.320074 

— 0,060109 + 0,028353 

+ 0.001210 

+ 0,000720 — 0.044537 

1830 

+ 0-323477 

—0.062456 +0.031106 

—0.007708 

—0.002107 —0,050635 

1840 

+ O.322334 

— 0,063180 +0,028641 

—0.000354 

— 0.001905 —0.048071 

f i860 

+ 0.325760 

— 0.059894 + 0,029350 

— 0.001917 

+ 0,00179s —0,047563 

1880 

'+ 0.333923 

— 0.061920 + 0.027656 

— 0,005507 

— 0.003805 — 0,051930 

1885 

+ 0.315720 

— 0,060258 + 0.0248 u 

+ 0.007906 

+ 0.012999 —0,049798 


years 

,(») „<2) 

-*» "a 

„( 3 ) 

•“0 

.0) 

■‘‘l 

B ?. 

. 13 ) „< 3 > 

. 4 , B j 

1538 

+ 0.02S063 —0.028945 

— 

-0.067942 - 

- 0.025209 

— 0,000341 —0.067016 

1572 

+ 0.035770 —0.01 1402 

— 

- 0.029237 + 0.00327 6 

— 0,010857 + O.OOO713 

1600 

+ 0021534 +0.011750 

— 

- 0.012933 + 0.016891 

+ O.OO1655 —0.0*1415 

1642 

+ 0.040967 —0.018373 

— 

-0. 1 14577 —0.032849 

— 0.008689 —0.008044 

1676 

+ 0.021315 —0.006077 

— 

-0058067 — 0.007136 

— 0.006700 — O.OI7533 

1700 

+ 0.004272 + 0.017030 

— 

- 0,003912 +0.007041 

+ 0,001426 —0.021382 

id, 

— 0.001686 +0.019182 

— 

— — 

— — 

1710 

+ 0.001774 + 0.016588 

— 

— 0008499 + 0.000764 

+ 0.000012 —O.OO95OO 

1720 

+ 0,006400 + 0.010162 

— 

—0.052379 +0.020753 

+ O.O08096 —O.O2544 5 

‘730 

+ 0.001351 + 0.012930 

— 

— 0,047195 +0.048S19 

+ O.OOO566 — 0,024583 

'744 

—0,002059 +0.018038 

— 

- O.O06576 + 0 , 00107 7 

+ 0.CO24S6 —O.OI7374 

1756 

— 0.004824 +0.014594 

— 

- 0,024272 — O.O08624 

+ O.Ol6llO —O.O573 95 

1760 

— 0.062791 — 0.015726 

— 

— — 

— _ 

1770 

— 0.0048 Z4 + 0.008602 

~, 

-0.01 1707 +0.027345 

— 0.006485 — 0,026892 

1780 

— 0.01 2609 + 0,0:6787 

— 0.016383 

- 0.019010 + 0.033642 

— 0.003230 — 0.020201 

1787 

— 0.006 133 + 0.007939 

— 

- 0.0 160 1 6 + 0,024783 

— O.O06418 —0.0277)8 

1800 

— 0.008586 +0.008619 

— 

- 0 026726 + 0.022926 

— 0,0038Z2 — 0.021547 

1820 

— 0,014513 —0.001898 

— 0,001792 

- 0.002657 + 0,027497 

— 0.005187 —O.OI3466 

1829 

— 0,012637 +0.001249 

— 0.01S763 

- 0,016446 + 0.029863 

— 0.tx62II — O.O3O7Z8 

1830 

— 0.013631 +0,000172 

— 0.006593 

- 0.016700 + 0,(^42956 

— O.X7955 — 0 , 0255/5 

1840 

— o.oii4s8 +0.000239 

— 0,016844 

-c.018081 +0.039628 

— O.Or6227 — O.O28549 

i860 

— 0.012562 — 0,002185 

— 0.026076 

- 0.010364 + 0.038904 

— 0,006550 — 0.03 1 145 

1880 

— 0,012433 —0.004284 

— 0.01079c 

- 0.003637 + 0.038543 

— 0 , 0 X 58 l — 0.032358 

1885 

— 0,012604 — 0.005667 

— 0.024363 

- 0.007383 + 0.039560 

— O.OX +43 —O.O27857 
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ORIGINAL PAG s' fs 
OF POOR QUALITY 


vWWi 

ft 


years 

,( 9 ) 

'*0 

d ? 

»S 38 

+ 0.006263 + 0,009728 

1572 

+ 0,004726 + 0,013087 

i6co 

+ 0,004264 + 0,001960 

164? 

+ 0.004333 + 0,019566 

1676 

+ 0,000957 + 0,014435 

1700 

— 0.001970 + 0,001662 

id, 

— — 

1710 

— 0.002055 + 0000248 

1720 

— 0.COS331 +0,006698 

, 73 Q 

— 0.005783 +0006873 

1744 

— 0.00655 1 + 0,002791 

756 

—0.009366 + 0,001333 

1760 

— — 

1770 

—0.008399 +0.001694 

1780 

— 0,004761 + 0,0079 1 8 

1787 

— 0.008470 +0.001649 

1800 

— 0.009338 +0.001366 

1820 

— 0.009141 +0.000736 

1829 

— 0007248 —0.001272 

1830 

— 0.006552 + 0,000488 

1840 

— 0,007175 - 

-0,001437 

i860 

— 0.006473 — 0.002960 

1880 

— o,cc68i2 —0.003239 

OO 

CO 

— 0.005492 - 

- 0.003270 


»«> ,(l) „U) 

■o 0 *1 1 « I 

— +0,204515—0.185687 

— + O.OI95I2 — 0,042398 

— + 0,035705 —0,041210 

— —0,245348 + 0.120C00 

— + 0.079424 — o, 1 1 8363 

— + 0,013895 — 0.101369 

— + 0.027658 — 0,005286 

— + 0,082152 — 0,097625 

— —0,054510 —0,045622 

— + 0,010207 — 0,016748 

— — 0.030747 —0.030085 

— + 0,003841 — 0,036963 

0.010504 . + 0,029979 — 0.059683 

— — 0,000563 —0,039763 

— + 0.038186 — 0064707 

0.002786 + 0.027844 — ^003377 

0.027377 —0,003197 “ 0,038405 

0,038035 + 0.022402 — 0,053291 

0,024255 —0.000065 — 0.033215 

0,023236 —0,005119 ->o, 0401 17 

0.018963 —0.014340 — 0.064692 

0,034395 —0.011877 — 0,030597 


years 

,< 4 ) n U> 

"■‘2 "2 

dr 

dr 

,u> 

Ai 

dr 

‘638 

+ 0.001029 + O.C93I58 

+ 0.012549 

+ 0,013520 

— 0.001437 

—0.002152 

‘> 7 2 

+ 0.026228 + 0.011724 

+ 0.003476 

— 0,006379 

— 0/200892 

— 0.000545 

1600 

+ O.O29696 —0.021 13I 

+ 0.001613 

+ 0.003107 

+ 0.002682 

— 0.000855 

1642 

+ O.O3I572 — 0.010402 

— o.coi 203 

+ 0.000303 

— 0.000905 

— 0,006608 

1676 

+ 0,016411 —O.O08635 

+ 0.005875 

— 0,002058 

— 0.001464 

-0.003795 

1700 

id. 

+ 0.020761 — O.OI5794 

+ 0.001076 

+ 0.000859 

+ 0.001026 

— 0.000026 

1710 

+ O.O24O/9 — O.OO2749 

+ 0,002019 

+ 0,0009! I 

+ 0.001692 + 0.001005 

1720 

+ O.O45469 —O.O36582 

+ 0.003260 

+ 0.004734 

+ 0.00068s 

+ 0,004087 

1730 

+ O.06C63O — O.066835 

+ 0.007790 

+ 0,007283 

+ 0.000818 

—0.001830 

‘744 

+ O.O26539 — O.OO7984 

+ 0.004447 

+ 0.001200 

+ 0.001 120 

+ 0.000448 

1756 

+ O.O3OI3I + O.COOI45 

+ 0.010959 

+ 0.022307 

— 0.000243 

+ 0,001258 

1760 

— — 

- 

— 

— 

— 

1770 

+ 0.024547 — 0,012719 

+ 0.006874 

+ 0.001956 

+ 0,001042 

+ 0.000290 

1780 

+ 0.0 1 0106 —0,000496 

+ 0.0020 1 6 

+ 0,005909 

+ 0.002973 

+ 0.002549 

1787 

+ 0.019176 —0.010822 

+ 0,007866 

+ 0.000485 

+ 0,001379 

+ 0,000171 

1800 

+ 0.016569 — 0,016654 

+ 0,003247 

— 0.001210 

+ 0.001691 

+ 0.000466 

1820 

+ 0,009462 — 0.035830 

+ 0,006754 

+ O.OO4I45 

+ 0.000908 + 0.001063 

1829 

+ 0,010851 —0,015592 

+ 0,003014 

+ 0,006708 

+ 0,001 109 + 0.000895 

1830 

+ 0,014876 —0.016000 

— 0.000062 

+ O.OC69O9 

+ 0.001109 

+ 0.001442 

1840 

+ 0.007561 —0,011004 

+ 0.003133 

+ O.O08165 

+ 0.000311 

+ 0.000149 

i860 

+ 0.017039 — 0,016568 

+ 0.004219 

+ 0,005815 

+ 0,000561 

+ 0.000336 

1 880 

+ 0.014439 —0.020615 

+ 0,0c 6755 

+ O.OIO284 

+ 0,001027 

— 0.000382 

OO 

CO 

Vi 

+ 0.007147 —0,019754 

+ 0,005121 

+ O.OC6842 

+ 0,000968 

— 0,000849 


/ a . 
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years weight «**> 


— 1 ' 

— ** ** 


A -A 

*538 

6 

— 

+ 0.071852 + '38*442 — 

y.103 ~ 

+ 0068652 + 94*.2i9 

+ 5 *. 83 i 

* 57 * 

6 

. — 

+ 0.084574 + 47 .687 — 1 

>4459 , . — 

+ 0.045493 + 180.417? 

— ■ 

1600 

8 

— 

+ 0.039203 + 33 .038 + 

3 -393 ' - 

+ 0.045418 +132.042 

- >5 -833 

1642 

6 

— 

+ 0.113758' + 54.937 - 

13.702 — 

+ 0041900.+ 132,179 

- 5.024 

1676 

6 

— 

+ 0.101950 + 47.706 — 

2 .382 — 

+ 0.033594 + 137.398 

— 1 .380 

1700 

8 

— 

+ 0.043474 + 41 .172 + 

6697 — 

+ 0.040695 + 146 .606 

~ 4-335 

id. 

i“/« 

+ 0.288251 

+ 0.094297 4 ' 39 - 4*9 - 

8450 —0,063404 

+ 0.087422 + 17,636? 

— 

1710 

7 

— 

+ 0.045608 + 44 .514 * 

4499 - 

+ 0.031053 + 152.230 

“ 7 -353 

17 20 

5 V« 

— 

+ 0.068197 + 53 .08a — 

2 - 9*3 - 

+ 0.046992 + 142 .531 

+ 4 - 95 * 

1730 

6 

. — 

+ 0.056559 + 53 .2*7 — 

2 -547 * - 

+ 0.049578 + 149429 

+ 0.661 

> 7+4 

7 

— 

+ 0.052866 + 51 .727 - 

1.17J - 

+ 0.048118- + 148.353 

+ 5.385 

1756 

57 * 

— 

+ 0.044486. + 55.943 - 

1 .669 — 

+ 0.094330 + 86,79 8? 

— 

1760 

s*/« 

+ 0.347964 

+ 0,068144 1-3.699? 

— + 0.024643 

+ 0.144615 +130,619? 

— 

1770 

9 

— 

+ 0.056660 + j6 .721 — 

0.845 ~ 

+ 0,040245 + 158.750 

+ 1 .764 

1780 

8 

+ 0.319216 

+ 0065151 + 43-507 + 

13.512 — o.oc6ioo 

+ 0.037464 + 1 66 .439 

- 3 3 *> 

1787 

9 

— 

+ 0.059247 + 58434 - 

’=.6.4 

+ 0.041823 + 157,895 

+ 7.050 

1800 

8 

— 

+ 0.062995 + 5 6 .869 - 

2,438 

+ 0.021927 + 152,279 

+ 16 .OJO 

1820 

21 

+ 0.333810 

+ 0.068341 +67.138 — 

5 143 +0.005092 

+ 0.030416 + 179,823 

- 6,278 

1829 

30 

+ 0.320074 

+ 0066460 + 64 .747 — 

2,123 + 0001210 

+ 0.044542 + 180,926 

- 5.03; 

1830 

21 

+ o. 3 2 3477 

+ 0069773 + 63.525 — 

• 0.787 —0007708 

+ 0,050678 + 177 ,617 

- 1.465 

1840 

"*7 

+ 0.322334 

+ 0.069368 + 65 .604 — 

1 .722 — 0.000354 

+ 0.048108 + 177.731 

+ 1 .027 

1860 

*7 

+ 0.325760 

+ O.C66697 + 63 .894 T 

2.276 — 0.001917 

+ 0 .C 47597 + ‘ 77-839 

+ 6.131 

jS3o 

30 

+ 0.333923 

+ O.067807 + 65 .947 - 

2 510 —0.00550 7 

+ 0.052069 + 175.809 

+ *3 -374 

1885 

75 

+ 0.315720 

+ O.065I67 + 67.618 - 

1 411 + 0.007906 

+ 0.0514 66 + 194 ,630 

— 4.144 


years j* *-A C ^ A- A « l .” fi* *- 

1 538 + 0.040316 —67’-944 — S’. 1 23 — . + c.c-2+70 +24^.644+ 8*441 +0.029373 + 89’,66; — 3’ 

1572 +0.037547—53.851—8.757 — + =.==420+276.393 — 9.061 + o.oic88o + 43 .121 + 41 . 

1600 +0.024531—30.690—20.834 — -=0:1273+322.559—47,612 +0.011534+94.117—9. 

1642 +0.044908—57.075 +22.177 — ■>•0.115193 +254.002+32.368 +0.011841 +68.603 +16. 

1676 +0,022139 — 52.957 + 31.518 — - ==r*j=3 + 262.994 + 32.623 +0.018770+79.543+ 5 

1700 +0.017557— 7-°4* — 4-898 — - =.^C5J + 339 .943 — 37 .799 +0.021430 +91.908— 6 

id. +0.019256 + 2.511 —14.450 — — — — — — — 

1710 +0.016683— 3.032--- 4.928 — - ccc?j33 + 84.863? — +0.009500 +90.000— 4, 

1720 +0.012009—16.104+12.082 — -=.=55=40 + 68.38 7? - +0.026702+98.825—13 

1730 +0.013002— 2.983 + *.920 — t 0.062847 + 44.031? — +0.024589 +90.659— 5 

1744 +0.018155 + 3.256 + 2.222 — - 0006663 +279.301 +34.810 +0.017546 +94.023— 8 

1756 +0.015372 + 9.146 + 1.083 — - =.=25-58 +289.560 +27.814 +0059613 +97.836— 12 

1760 +0.064730 +52.039? — — — — — . — — — 

1770 +0.009862 +14,642 + 1.129 “ -=.=25+45 + 336.822 — 15 .640 +0.027663 +§3.221 + 2 

1780 +6.014320 +30.854 — 11.125 —0.016383 -=.::?64X +33 0.531 — 6,629 +0.020457 +85458 + o 

1787 +0.010032 +18.843 + 3-^57 — -==55=8 +327.127— 1.322 + 0.028491 +83490 + 2 

1800 -0.012165 +22.446 + 5.200 — - =.=55212 + 316.193 +13.148 +0.021883 +84.975 + o 

1820 +0.014637 +48.725 — 13.162 —0,001792 +0:2*625 +354.481 —19.701 +0.016312 +80.729 + 5 

18Z9 +0.012699 +42.177 — 3.051 —0.018763 -0:52=85 +331.152 + 6.076 +0.031349 + 84.2S6 + 1 

1830 +0.013632 +44.638— 5.116 —0.006593 + =.=6087 + 338 .755 — 1.255 +0.026783 +81.361 + 4 

1840 +0.011461 +44402— 0,922 —0.016844 +0.225558 +335475+ 4.745 +0.029221 +83.848 + 2 
i860 +0012751 + 49-933 + 1484. —0.026076 — : 22=61 +345.083 ; + 0.576 +0.031827 +84.061 + 1 

1880 +0.013150 + 54.506 + 4.808 —6.010790 - :.:;.**13 + 354.609 — 3.510 +0.032362 +89,485 — 3 

1885 -0.013819 +57.104 + 4.189 —0.024363 - : '-:=3 +349.429 * 3.029 -‘-0.027860 +89.544— 3 


ORIGINAL PAGE & 

OF POOR QUALITY 


years 

(3) 

«3 

PT 


~Po 

‘538 

+ 0,012658 

- 9 \885 

— 

3’.202 

1572 

4- O.OI3914 

— 6,618 

— 

t .461 

1600 

4* 0,004691 

— 21 ,768 

+ 17 813 

1642 

4* 0.020040 

■■ * 4 *164 “t* 

6*395 

1676 

4* 0.014467 

— 1 .264 

+ 

8,503 

1700 

4- 0.002577 

4- 16.616 

— 

5,842 

id. 

— 

— 


— 

1710 

+ 0.002070 

4- 27 .706 

— 

*5 '459 

1720 

4- 0.008561 4- 12 .839 4* 

0.881 

•730 

4- O.OC8982 4* 13 .350 4- 

3,834 

1744 

4* O.CO7121 

4- 22 .308 

— 

5.053 

1756 

4-0,009163 4*27,212 

— 

8.190 

1760 

— 

— 


— 

1770 

4* 0,008568 

4- 26.199 

— 

5,115 

1780 

4* 0,009241 

+ 10,371 

4* 12,186 

1787 

4- 0.008629 + 26 -328 

_ 

2.740 

1800 

4* 0.009140 4- 27,135 

— 

1 .632 

1820 

4- 0.011365 4- 17.847 

4- 10,602 

1829 

4- 0.007359 + 33-3> 8 

— 

3-544 

1830 

4-0.008170 4- 17 774 + 12.148 

1840 

4- 0,007318 

+ 33 775 

— 

2,380 

1860 

4- 0.007117 4- 38.191 

— 

3 .85» 

1880 

4* O.O07543 

4-38 .477 

•— 

i ,191 

1885 

4- 0,006392 

+ 40,257 

— • 

2.234 


rt (4) (i) ^»4) a a 

«0 «t t*l PC~~ Pq 

— > + 0.276240 +■ zzfj 6 z — 3’, 205 

— 4 - 0.046672 + 204,712 + 15 .072 

— + 0-254525 + 220 .906 — 5 ,053 

— 4-0.273121 + 26.064? — 

— 4*0.142539 4-213.863 — 8,680 

— 4*0,102316 4 * 187,805 4- 14.009 

— 4- 0.028159 4* 259.180 — 58770 

— 4*0,127591 4- 220 ,c8i — 21 .075 

— 4-0.071081 4* 129.929 4* 67,673 

— 4- 0,019613 4- 21 1 .360 — 15 723 

— +0,043017 4 - 134.377 + 59-575 

— 4* 0,037162 4*185.932 4- 6,055 

— 0,010504 4 * 0,066790 4 * 2 c 6 ,671 — 16 ,088 

— 4*0.639766 4* 180.810 4* 8,790 

— 4 - 0,075134 4-210 .548 — 22 773 

— 0,002786 4* 0,028048 4*^263 .085 — 78 ,1 18 

— 0.027377 4-0.038536 4-175-241 + 8.463 

— 0038035 4-0,057809 4-202,801 —19,238 

— 0,024255 4*0,033215 4*179,888 + 2,271 

— 0,023236 4-0,040442 4 * 172728 4 - 6,624 

— 0,018963 + 0,066261 4 * 167.502 4 * 9 ,042 

— 0,034395 4-0,032819 4* 158.784 4- 17,058 


/21 


years 



(*% 

P * —A 

4 * 


Pc ~~ Pv 

«V' ft ' Pc-Po 

' ‘538 

4 * 0,993162 

4 * 

i 79’.683 ? 

— 

4- 0,018444 4- 345’.7*3 

— 

Z 3 *. 2 3 S 

4- 0,002587 4 - 36’. 567 4 * 13*212 

1572 

4 - 0,028729 

4 - I47.042 

4- 

6’.375 

4- 0,007264 + 309 ,529 

+ 1 5 .982 

4- 0,001045 4- 30.356 4* 

21 .435 

1600 

4 * 0,036447 4 - 117,282 

+ 18,426 

4 - 0,003501 

+ 3 jo .855 

— 

22 ,847 

+ 0x02814 4 - 63 .077 ? 

— 

1642 

4 * 0,033148 

4 - 

126.128 

4 - 

4-594 

■ 4 * 0.001241 

4- 25 ,288 ? 

— 

4- 0.006670 4- 43 .051 4* 12 ,884 

1676 

+ 0 , 0185*14 

121 ,124 

4 - 

5.561 

4 - 0.006225 

+ 323 .565 + 

11 .221 

4 * 0.004068 4 - 39.726 4 - 18 .222 

1700 

id. 

+ 0,026086 

+ 

116.369 4* 

7.467 

4 * 0,001377 

4 * 342 .867 

— 

5,940 

4-0,001026 4-67.137 — 

7.768 

1710 

+ 0,024236 

+ 13 1 -743 

— 

9,094 

4 - O.OO2215 

+ 338 .095 

— 

OX76 

4 - 0,001968 4 - 75.177 — 

15 .216 

J720 

+ 0.058356 4 - 115.590 

4- 

5.872 

4- 0,005748 4* 348 ,482 

— 

9 77 ‘ 

4-O.OO4143 4-87.635? 

— 

> 73 ° 

4- 0.090238 

4 * 

Mi .ioS 4- 

9,168 

4 - 0,010664 + 344 .358 

— 

4*755 

4 - 0.002004 4 - 51 .021 4 - 

10.124 

*744 

4 - O.027713 

4- 126.628 

-- 

8,016 

4- 0.004606 4- 335 .034 

4- 

S.817 

4- 0,00 1 zc6 4- 72 .952 — 

10 ,978 

1756 

1760 

1770 

4 - O.030131 

+ 

179,862 ? 

— 

4 * 0,024853 4- 3 >I .279 

— 

9.357 

4-0.001281 4-92.733? 

— 

4 - 0.027646 

u. 

12! .354 

— 

5.828 

4 * 0,007147 

+ 335 ’ 2 95 

— 

7.875 

4- 0.00108 1 4- 71 .378 — 

7 .865 

1780 

4 * 0.0101 18 

T 

'33 -595 

— 

19,257 

4 * 0,00624! 

+ 353 720 

— 

9.638 

4- 0,003023 4- 70,115 — 

6 ,010 

1787 

4- 0,022019 

+ 

120.280 


6.773 

+ 0,007881 +331.176 4-13.511 

4- 0,001787 4- 69.267 — 

4.748 

1800 

4- O.023492 

+ 

H 2-573 

— 

0.609 

4* 0,003463 

4- 323.187 4- 

22 .659 

+ 0x01754 + 7 1 .352 — 

6x63 

1820 

4* 0,037058 

*r 

97 .396 

4 * 

12,194 

4- 0x04210 

+ 356 .644 

— 

9,014 

4- 0 COI 398 4- 79.873 — 

13 ,400 

1829 

+ 0.018996 *f 107.418 

+ 

i .103 

+ P .007354 

+ 35 ' -935 

— 

3.502 

4-0.001425 4-77.226 — 

10,220 

1830 

+ 0.021847 4* 

11 1.457 

— 

3.055 

+ 0.006909 4* 0.171 

— 

1 1 ,649 

4 - 0.001819 + 80,609 — 

>3 .544 

1840 

+ 0.013351 

*r 

107,24; 

— 

0.032 

+ 0x08745 

+ 353 .°°3 

— 

3.589 

+ 0.000345 + 73 .900 — 

6.243 

i860 

+ 0.023766 

- 

112.901 

— 

8.c6o 

4- 0,007184 

4- 348 .012 

4 * 

3 .'85 

+ 0x00654 4- 75 .230 — 

6.389 

1880 

4* 0,025168 

+ 107.504 

— 

5. 0 38 

+ 0,012304 

+ 348 .900 + 

4 .08 ! 

+ 0.001096 4- 6z .399 + 

7 .626 

OO 

OC 

+ 0,021007 

+ 

99*945 

4 * 

1 .928 

+ 0.008546 

+ 347.729 

*f 

5.698 

+ 0.001288 + 57 ,187 + 

'3.134 
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below. The y^ values correspond to the year 18 00. 


«o = + 0.322946 

K t 11 as 4. O.06705J 

n « a + 0,00K06- 
„ia) , : 

«1 =*+0,047815 
n t - + 0,013244 

Jt) 

a., 0,0,75:27 

, «t ~ + 0.030018 

(3) 

, ** + °' 028 373 
n i — + 0,007474 

=~ O.O26447 
a i - + O.040680 

* < t ) 

rt * , — + 0,022245 
"s = + 0,008200 

„U> 

a * ~ + o.ooj 1 55 


ji) 

~ 59>P7 


^ = < 68.553 

= 2-646 


,/tSt 

y 1 = 
y't * 


5=9-344 

85,688 
=2 J03 


?U • 

/ 2 = 1 1 1 564 
/a** ~ 34 ; 846 
//' = 65 289 


»>l° a* +• Cj’,11438 

«», 1 si 4- 0,26062 
«<**’ a + 0,39585 

„,<3> 

*» + 0 , 27(97 
w* =5+0 .00540 
Wj’ 1 as 4. 0,14729 

,„U> 

u ““-°' , 4°39 
w a « — 0,11872 
wi 3 =4-0 ,08919 
Wt 4 SI 4. 0,05920 


By calculating these numbers we have given an un- 
equal weight to the different years and we estimated 
this weight on the basis of the number of parallels 
entering in the condition equations. The weight of the 
incomplete observations of longitude decreased by a 1 
ratio of 3/4, based on the number of meridians shown on 
the charts. I allowed myself to exlude a small number of 
observations which exhibited extraordinary discrepancies. 

(/Is ) 

^ n s ^ us eliminated are followed by a question mark. 
The a til coefficients were calculated on the basis of the 
most accurate data of the present century, where the in- 
tensity of the magnetic force is known. Furthermore, I 
found that even if i included ancient observations, the 
values adopted do not change appreciably. 


The errors resulting from these formulas in the 

(% ) ‘ \ 

$ n angles are shown -hext to thes& in the 8-6 columns. 

' c o 

If we Inspect these numbers, which are often quite large, 

we cannot say that our corrections gave completely satis- 
factory results. Yet, most of the functions move in 

YL 


a positive direction, like the lower order terms. Only 
two of them [the Y ^ and Y functions] exhibit a nega- 
tive direction. Furthermore, 1 still have some reserva- 
tions about this, because these functions are certainly 
very incorrectly determined, as may be concluded from 
the enormous values acquired by the $ c — B 0 deviations. 

In summary, we arrive at this somewhat discouraging 
result , that from the seven new functions included in 
our second approximation, only one or two seem to be de- 
termined with sufficient accuracy: these are the third 
order functions [ Y and Y^]. 

Without any hesitation, we ascribe this result to 
the Imprecisions in Hansteen's charts. As a matter of 
fact, we have based all of our calculations on the carto- 
graphic representations, which we sketched by hand, ac- 
cording to documents which are often highly precarious, 
and by allowing ourselves to be guided by a sense of 
continuity. This was obviously a very unlogical way to 
correct the observation errors, but it greatly simplified 
the calculations. 

Now that we have approximated values of the unknowns, 
we can advantageously substitute this for a more rigorous 
method of calculation. This work must be preceded by a 
very thorough comparison of the theory with the observa- 
tions: each observation of the declination, inclination 
or intensity, made at any period of time, -will result in 
a difference; observation minus theory, and the set of 
corrections thus found will make it possible to obtain 
the most probable values of the unknowns and to thus 
reach highest level of perfection of which the theory 
seems capable. 
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Although we must not consider the numerical part of 
our work to be completed, the results obtained already 
make it possible to formulate certain statements concern- 
ing the general form of the secular inequalities. 

Our calculations bring us to represent the function 
of the magnetic forces at' a' single given place with a 
series of terms having the expression: 

JP® cos (mbt + nyf + nm® t), 


where the term t stands for the cosine sign and is the 
only variable element. A similar form is convenient for 
representing the three components of force, except for 
the component perpendicular to the meridian, the cosine 
must be replaced by a sine (no. 6). 

The expression for the horizontal component perpen- 
dicular to the meridian does not contain a constant term: 
for a long enough time factor, its mean value is zero, 
and the value of the component oscillates about the zero 
value. This must also be true for the magnetic declina- 
tion, at least approximatively . 

The other two components of force, which are within /2 4 

the meridian plane, contain in addition to the variable 
part, a constant part about which the horizontal and 
vertical intensities oscillate. 

Besides these constant parts, each of the components 
contains a series of periodical terms with a more or less 
long cycle. 
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If we continue to let be the average yearly longi- 

n (i) 

tudinal displacement of argument 0 ‘ , the revolution 

/ • j o n 

time of the term ' will be ~~u) and this term will 

n 

generate a corresponding cycle term 2 ir , . , , in the expres- 

nm n 

sions of the components of force. 


I find in particular, for the main terms for which 

calculations already made make it: possible to determine 

(i J 

the yearly motion m with some accuracy, the following 

Tv 

cycle values: 


Term 


,-U! 
1 i 
,.( 3 ) 
1 3 


Cycle 
x, r, z 

., 4 -Yrs Yrs 
15S1 1381 

002 4.C.1 

r < tJ | 

* 3 — *324 

2+12 815 


For the elements provided directly by observation, 
the theory does not indicate any simple expression, and 
it is easily conceivable that there are no determined 
cycles for different places, as we often thought we 
noticed . 

The magnetic axis of the earth is the direction of 
the main magnetic moment of the earth. With this, we 
find that the earth's magnetic axis descibes a circular 
cone about the rotating axis, with an angular distance 
of 23° 28' in 3147 years. According to the results of our 
calculations, the magnitude of this moment is invariable 
and equal to JR 3 . Va ^ 2 4 , , or in an absolute value, to 

the number 852640 followed by 20 figures, if centimeters 
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grams, seconds for average time are adopted for units. 

The displacements of the earth's magnetic poles 
are not as simple. These points are characterized each 
instant by the condition X « o, Y « o, where X and Y 
are in the form mentioned earlier. These are two tran- 
scendental equations in y and after eliminating time, 
they give an equation between y and which defines a 
cone whose outline with the sphere is the trajectory 
described by the magnetic poles over the earth's surface. 

If we limit ourselves to the first order terms only, 
the magnetic poles would become indistinguishable from 
the points where the magnetic axis passes through the 
terrestrial surface and would move with them. Moreover, 
the higher order terms essentially change the nature of 
the trajectory. If we again adopt the uniform motion for 
the first approximation, we obtain equations with finite 
differences which we satisfy by introducing the periodi- 
cal terms in y and w. We may thus account for the gener- 
ation of these yawed curves which various authors have ob- 
tained using a completely different method ( Astronomlsche 
Nachrichten , No.s 3254 and 3299.),* 

The earth's magnetic equator satisfies for each of 
these points, the conditions Z = o, where Z is a known 
function of y and This is a complicated equation which 
would seem embarrassing to discuss. The points of inter- 
section with the terrestrial equator (y = o) , are easier 
to determine. The longitudes of the ascending and descend- 
ing nodes of the magnetic equator on the terrestrial equa- 
tor are given by the following equation: 


*We read in most treatises that the earth's magnetic poles 
describe a complete circumference in about 600 years i this 
result is not accurate: we have found that it is the sec- 

ond and third order terms which make it seem like a cycle 
of about 600 years. 
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o = ~ <4 W + 2f< [ l) cos (nr + /j U 4. 

+ 3«? cos (star + zyf + aw*/) 

— cos (nr + yW + 

+ 4«f c °s (3W + 3y? + 3 ,n ? t ) 

+ f< 

— I «”’ COS (2RT + + 2U\^H) 

+ S 0 ?’ cos (4W + 4/" + 4 wjty). 

This is an equation of the eighth degree in cos ttr 
with real coefficients. For each determined value of t, 
the equation with eight roots including one or several 
pairs may be real. 
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We had nothing to show in these pages for a theory 
with a numerical perfection, which can be achieved only 
by operating step by step. Our objective was more modest. 

This was to show a method by which we can hope to 
penetrate into a difficult subject and to shbw the re- 
sults which we may hopefully obtain through a systematic 
application of calculus: this is the task which we were 
allowed to undertake and which we tried to fulfill. 

However, in the results already obtained, we will 
find solid ground upon which we may base all of our 
subsequent investigations and will enable us to calcu- 
late with the precision we want the higher order terms 
which we missed during the first test. 

A comparison of the theory with the observations will 
give us a better idea of the direction to take in such a 


% 


* 


revised calculation. It is for this purpose only that we 
want to assess the still imperfect results of a comparison 
of the theory with observations. 

For our comparison term, we have selected the declina- 
tion and the inclination for the year 1700. For this date, 
we have calculated the magnetic force and its direction 
based on ;the numerical constants obtained in the preceding 
pages. 


We will succinctly indicate the steps taken in the 
calculation of the magnetic force and its direction. 


/VI (A) 

The modules a of the functions f and their argu- 

* a ) n n 

ments 0^ are first calculated according to the formulas 
of no. 9, by using the numerical values given in this chap- 
ter. 


As these transcendental quantities are assumed to be 

U) J-i) 


known, we will easily pass to the constants A' 1 " and B 

n n 

using the formulas of no. 5, after which we may easily 

obtain the coefficients L . L’j M t M’j N. N'j which 

n n n n n n 

appear in the expansions of the forces as a function of 
the lorigitude multiples. . These coefficients, and the 
components of force X» Y and Z themselves, are given by 
the formulas of no. 6. 


/26 


Now all we have to do is calculate the elements we 
are lookin for; S, i and £ which are obtained using a 
well-known system described in no. 1. 


The tables below contain the intensity and the direc- 
tion of the earth’s magnetic attraction in the year 1600, 
calculated on the basis of a provisional equation for the 
points of intersection of thirteen parallels egually spaaed 
with twelve meridians. 
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Total force 


0 

or ho 1 

3 °’ 

60’ 

9 °’ 

I20 1 

• 5 ° 

180* 

zio* 

240’ 

270’ 

300’ 

33 °’ 

0 

0.5919 

0,5919 

0.5919 

0.5919 

O.59I9 

0-5919 

0.5919 

o. 59»9 

0.5919 

0,5919 

0,5919 

0,5919 

'5 

0,5728 

0.5749 

0.5822 

0,5904 

0.5949 

0.5938 

0,5921 

0-5955 

0.6:05 

0,6001 

0.5909 

0.5789 

3 ° 

0,5464 

0.5642 

0,5729 

0,5665 

0,5549 

0-5459 

0,5565 

0.5902 

0.6227 

0.6203 

0.5828 

0.5480 

45 

0,5098 

O.S 259 

0.5223 

0.5020 

0,4813 

04673 

0.4819 

0.5433 

0,6149 

0,6284 

0.5743 

0.5(69 

6o 

0.4575 

0,4581 

0.4344 

0,4063 

0.3912 

0.3857 

0.4023 

0.4686 

0.5604 

0.5893 

0.5356 

0,472 ( 

75 

0,3820 

0.3702 

0-3344 

0,3406 

O.3363 

0.3387 

o .3474 

0,3908 

0,4761 

o. 5”4 

0.4679 

0.4071 

90 

0.3(72 

0.3034 

0.3156 

0.3287 

O.34O9 

0.3589 

0.3585 

0,3725 

0.40S5 

0.4366 

0,4071 

0.3516 

105 

0.2907 

0,2784 

0,3109 

0 . 35 H 

O.3852 

0.41:3 

0.4009 

0.3770 

0.401 1 

0.4328 

0.4091 

0,3429 

120 

0,3172 

0.3039 

0.3503 

0.4178 

O.474I 

0.50-3 

0,4919 

0.4597 

0.4755 

0,5078 

0,4799 

0-3933 

*35 

0.3853 

0.3655 

0.4097 

0.4835 

O.5525 

0,5894 

0.5796 

0.5561 

0.5651 

0.5833 

0.5501 

0.4642 

150 

0,4763 

0,4600 

0,4936 

0-5553 

0.6153 

0,6494 

0.6507 

0.6402 

0.6397 

0.6367 

0.6016 

0-5357 

165 

0.5756 

0.5685 

0,5846 

0,6159 

0,6489 

0.6717 

0,6804 

0.6793 

0,6729 

0.6595 

0.6352 

0.6027 

180 

0,6526 

0,6526 

0,6526 

0.6526 

O.6526 

0,6526 

0.6526 

0,6526 

0.6526 

0.6526 

0.6526 

0,6526 


Inclination 


0 

rfso’ 

30 * 

60’ 

90- 

120’ 

150’ 

180’ 

210* 

240’ 

270’ 

300° 

33 °’ 

0’ 

+ 86' 13’ 

4 - 86’ 13' 

+ 86" 13' 

+ 86’i3' 

+ 86-13' 

+ 86’ i 3’ 

- 86-13' 

+ 86’13' 

+ 86’iy 

+ 86’ 1 3' 

+ 86’ 1 3' 

+ 86*13; 

15 

4 - 84 15 

+ 8344 

4 - 82 26 

+ 8018 

+ 78 29 

+ /8 3 

-79 3 

+ 8125 

+ 8515 

+ 8927 

+ 87 22 

+ 8s 19' 

30 

4- 78 36 

+ 77 4 ' 

+ 7440 

4 - 71 2 

+ 6751 

4-6623 

-5739 

+ 73 1 

+ 8034 

4 - 85 24 

+ 82 16 

+ 7924 

45 

4 - 7032 

4- 68 1 1 

4- 62 29 

+ 57 4 

+ 53=7 

+ 5022 

“ 51 jo 

4-6056 

+ 71 58 

+ 7719 

+ 7447 

+ 7139 

60 

+ 5822 

+ 54 ' 

4 - 4440 

+ 37 3 

+ 32 54 

+ 2731 

-29 5 

+ 4+4 

+ 59 *3 

+ 6+53 

+ 62 55 

+ S 95 « 

75 

t 42 35 

+ 34 46 

4 - 1836 

+ 83 

4 * 436 

+ 138 

- 046 

4 - 21 47 

+ 3945 

4 - 4421 

+ 43 9 

4 - 42 42 

90 

4 * 1948 

4-916 

— ll 25 

— 2223 

-2432 

“29 39 

—•2546 

- 5 36 

4 - 11 52 

4 - 1424 

+ 13 47 

4 - 17 2 

1=5 

- 423 

“ 134 * 

“33 11 

-4252 

-44 52 

— 48 17 

— *5 4 

“29 5 

— 1813 

— 18 50 

— 1829 

— II 13 

120 

— 31 22 

— 3623 

-4918 

“5738 

— 6036 

— 62 4.3 

“59 23 

“So 7 

“45 5 * 

“47 54 

-4647 

-38 56 

'35 

“47 59 

-4927 

“57 55 

-65 39 

-7027 

-72 32 

“^24 

— 64 11 

-6317 

— 65 20 

-63 4 

“55 6 

150 

-6035 

-6019 

-64 53 

-71 17 

-7720 

-8055 

--313 

-76 8 

—7536 

“75 53 

-7225 

-6550 

165 

— 71 12 

-7047 

“ 7 * 5 * 

1 

NJ 

Ov 

«*- 

— 81 25 

-86 7 

“ 5 * 5 ! 

— 8625 

— 84 22 

— 81 48 

— 78 2 

“74 0 

180 

-8l (8 

-81 (S 

-81 18 

-81 18 

-81 18 

-Si iS 

— 1 1 rS 

— 81 18 

-81 18 

-81 18 

-81 18 

-81 18 
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Declination 


Q 

ter » o’ 

3 °’ 

60 5 

90 ’ 


120’ 


150* 1 So” 

210'' 

240' 

270’ 

3 CO 1 


330’ 

O’ 

+ »I 5 , I 3 ' 

+ w 

+ 55 ** 3 ' 

+ 2,-13' 

— 

4 * 47 * 

— 

3 ^ 47 ' — 64 , 47 ' 

“• 94 * 47 ’ 

-I 24 > 47 ' 

-1 54 * 47 ' 

+ » 7 S*' 3 ’ 

+ 45 ’* 3 ' 

'5 

+ 27 9 

+ 2537 

+ 2538 

+ *8 33 

+ 

7 * 

— 

8 13 - 2743 

- 48 39 

- 6852 

~ 044 

+ 6648 

+ 41 20 

30 

+ 

1047 

+ 13 1 

+ 1647 

+ «4 7 

+ 

732 

— 

2 +f - *747 

- 33 47 

— 41 10 

— 2 12 

+ 3055 


21 1 

If 


711 

+ ii 22 

+ *4 9 

+ 1053 

+ 

S 3 * 

— 

038 — 1352 

— 2651 

- 27 5 

- 2 5 

+ 1730 

+ 

* 3*9 


+ 

6 3 

+ 11 40 

+ > 33 * 

+ 9 5 

+ 

534 

+ 

029 — 11 S3 

— 2238 

— 1921 

— 228 

+ 8 58 

+ 

743 

? 7 S 

S 

528 

+ 1249 

+ '534 

+ 811 

+ 

5 5 

+ 

1 8 — 11 32 

— 20 3 

- 1346 

- 217 

+ 3 20 

+ 

3*4 

\go 


4 * 

+■ »4 3 

+ *552 

+ 9 8 

+ 

5 4 ' 

+ 

123 - 1*43 

- 1729 

~ 930 

- 158 

— 041 

— 

* «3 

°S 

+ 

1 23 

+ *4 55 

+ 19 18 

+ 1240 

+ 

749 

+ 

1 18 — 10 14 

~ '443 

- 5 56 

— 1 46 

- 430 

— 

6 Ik 

12° 


141 

+ 1417 

+ 23 11 

+ ' 9*5 

+ 

1220 

+ 

r 42 — 12 9 

- 1157 

— 326 

— 23* 

- 93 ' 

— 

12 15 

hs 

— 

645 

+ 12 0 

+ 2 55 * 

+ 27 59 

+ 

2030 

+ 

3 54 ~ "37 

— 1031 

- 328 

- 6 5 

— 17 0 

— 

1840 

ISO 

— 

940 

+ 10 8 

+ 27 2 7 

+ 37*6 + 35 4 8 

+ 

>3 53 - a 4 

- '257^ 

95 z 

— 17 26 • 

— 28 7 

— 

2442 

165 

— 

11 6 

+ 1023 

+ 3134 

+ 5048 

+ 6639 

+ 

75 25 - 5 'z 

~ 4547 

- 4851 

— 5045 

- 4444 

— 

30 26 

180- 

i 


11 1 

+ 18 59 

+ 48 59 

+ 78 59 

+ 10859 

+ 13859 +16859 

—161 1 

—131 1 

— 101 1 

.- 7 ' ' 

— 

41 1 
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Comparison of the theory with observations 


0 

VT =* o’ 


3 °* 

60’ 

90’ 

120’ 

'S'? 5 


I8rf 

2 1 o’ 

240- 

2 7° : 

300’ 


33 0? 

30 ' 

+ I 9 ’.z 8 

+ 1 2’, 42 

— o', or 


+ 7*54 

+ if ,27 

+ 22’, 22 

-26 1 , 79 

— ioi’,i6 

— 39’.20 

+ 9’.32 

+ if. 72 

45 

+ 13 .99 

+ 

9.16 

+ 1 .25 

— 1 ,61 

+ 3.32 

+ 12.37 

+ 

16.I3 

— 21 .85 

— 47,08 

-24,69 

+ 2 ,50 

+ 

12.61 

60 

+ 875 

+ 

6.46 

-1.67 

— 1 .11 

+ 5,96 

+ 11 .48 

+ 

8,12 

~ 12.63 

- 14-65 

- 16 ,86 

+ 0,36 

+ 

9-53 

75 

+ 6,67 

+ 

5 - 3 ' 

“ 3-03 

— 1 .82 

+ 7 .08 

+ 11 ,14 

+ 

5-47 

“ 9 '°S 

“ 9 7 7 

— 8.9S 

+ 0,14 

+ 

6.93 

90 

+ 6.63 

+ 

6 ,25 

- 5.44 

— 2 ,8 7 

+ 7 . 09 

+ 10 .39 

— 

I .72 

— 11 48 

— 6.50 

— 3.96 

— 0,89 

+ 

.7.18 

105 

+ 8.09 

+ 

8.£3 

“ 3 - 9 ° 

-3.04 

+ 6,12 

+ 7 79 

— 

4 >z 3 

-14.72 

- 6.43 

— 3.60 

— 0,51 

+ 

7.12 

120 

+ 6.71 

+ 

9,28 

— 0.42 

— 1 .76 

+ 5.14 

+ 5.69 

— 

5.65 

-12.95 

- 7.44 

- 3. S' 

+ 0.49 

+ 

6.94 

*35 

+ 6.95 

+ 

IO . 20 ’ 

+ 4.26 

+ 0,09 

+ 7.10 

+ 6.90 

— 

3.6' 

— 7 . 5 2 

- 3 48 

— 2.38 

-3,40 

+ 

6.64 


The last table contains the results of the compari- 
son of the theory with observations. In it we show the 
differences between the calculated and observed declin- 
ations, in direction & c <$ q . The pattern of the dif- 

ferences shows that they are essentially due to the terms 
multiplied by the arguments of the expression 2 vs + 2 yf + 2 mft. 


ORIGINAL’ PAGE F3 
OF POOR QUALITY 

A quick calculation has given the following results for 
the arguments of the terras dependent upon o' and aw, in 
the third approximation: 


A l> - 37’M9, 
iiT - ni .80 5/ 
Pi* = 266.631, 
/i U = >3 798 ; 


A" « 3^893, 

-a 290,256, 
^=283.765,. 


from. -Which we would conclude that the motion of all fourth 
order arguments is fairly rapid, in A positive direction . 

There are a certain number of inclination observations /28 
at the turn of the XVII th centry, made by scientists in 
Europe and by English voyagers looking for north-west and 
north-east passages. We are reporting these observations 
here, although may be considered only very rough figures. 

We have made important corrections to the inclinations 
of Baffin, one due to a collimation error, and the other 
to an incorrectly balanced suspension. The rest of the 
observations were left the same. 

We may find it of some interest to compare them with 
the theory, and we can do this appoximately by using our 
inclination table. We will find that the theory represents 
the observations to the nearest 1 or 2°, and as the incli- 
nation has diminished by about 6° over the years, we see 
that the theory already accounts for most of the variation. 

To have a better appreciation of the value of this result, 
let us remember that the oldest inclination observations, 
which served to define the arbitraries of the theory, were 
borrowed from Hansteen’s chart for 1780. 


41 


ORIGINAL PAGE Is 
0F POOR QUALITY 


Observations Of The Magnetic Inclination 
In The XVI th And XVIIth Centuries 


a, 

• Licux. 



b 

C 

d 

e • 

f 



Latitude 

cp—O 

Long; rsie 
17 

Dates. 

, 23 

Incllnaison 

i 

Autoritds. 

1 Hornsound .... 

• 

. , 

+ ‘fi’-ss' 

I5.5- 

1613 jum 14 

4 86’. 0' 

Baffin. 

2 Mer glaciate . . . 


» * 

+ 75 - zz 

41 , Si- 

t6o8 juin 19 

4 89 .30 

Hudson. 

» i . « 


, , 

+ 74-35 

44 

1608 juin 22 

4 85 , 0 

Hudson. 

V f • • 


, , 

+ 74 - z 3 

26. 

t 6 o 3 juin 7 

t 86 • 0 

Hudson. 

3 Swnrle-Cliff (M'ouv. 

Zcm e ) 

-f ~z ,iz 

51. : 

16=8 juin 29 

T 84. C 

Hudson, 

4 Cap-Nord .... 


♦ ♦ 

+ 71 , o 

25 . ;r 

1608 j u= 2 4 3 

+ 84 ,30 

Hudson. 

5 C6le de la Norv&ge 


• * 

4 69.40 

17 -45 0 

t6c8 nl’at 26 

+ 84 ; 0 

Hudson. 




-r 67. O 

« 3 - ? 

1608 mat 22 

+ 82 . O 

Hudson. 

6 gale de Baffin . . 


• # 

+ < 5 j .45 

3 S 3 • c 

1613 mai 23" 

+ 82 . C 

Baffin. 

7 C6le de la Norvige 


* t 

4 64.52 

to 

1608 

+ 81 . 0 

Hudson. 

8|Frobi3her Strait . . 


i « 

4 61 .40 

295. : 

1602 jdiT.L .7 

r _ 24 

1607 mat 3 

+ 83 .30 

Weymouth, 

9 Pres Shetland . . 



t- 61 .11 

35 “ 4 = 

4 79 - 0 

Hudson. 

lOlfranibourg .... 



+ 55 -54 

1242 

j 583-1585 

+ 73 - 2 5 

Tycho- Brahd. 

llLondres . . . . . 



+ 5 ' - 3 » 

359 - 5 ' 

1576 

4 7 « -50 

R. Norman. 

> , , . . . 



t 


1600 

4 72 . 0 

Gilbert. 

1 « « » t » 



V 


1613 

4 72 .30 

Ridley, 

t , , , , . 



* 


1676 

+ 73 - 3 ° 

Bond. 

12 Rouen 


• • 

+ 49 .26 

1 . f 

' 639 > 

4 7 z - 0 

P. Grandatni. 

13 Paris 


♦ • 

+ 4S .50 

2 2: 

l6 “‘ 26 

+ 75 • 0 

Richer. 

14 Kauffbeuren . . . 


• » 

+ 47 -53 

12 52 

1607 nov. 28 

- 70. 0 

Kdplcr. 

15 Tours 


t * 

+ 47 .24 

55 ? i* 

1659 :• 

4 70 , 0 

P. Grandami. 

16 Ferrare 


» » 

+ 44 - 5 * 

it 5- 

1629 

4 62. 0 

Cabreus. 

17 Rome 


• • 

4 * 4 1 -54 

12 2- 

■1639 '■ 

4 65 .50 

P. Kircher. 

18 Maples . . . . , . 



4 - 40 .52 

U 12 

1659? 

+ 59 1 5 

P. Kircher. 

19 Lisbonnc 



+ 40 . 0 

552. 1' 

1639 ? 

4 62 ,45 

P. Kircher. 

20 Messine 



T 5S.II 

15 ■: 

1639- 

4 59.15 

P. Kircher. 

21 Palermo 



+ 38. a 

15 -• 

1639 

4 59.15 

P, Kircher. 

22^‘e 


• • 

+ 35 -54 

15 

1659; 

4 59 .i; 

P, Kircher. 


Key: a -Place:; b-Latitude; c-Longitude; 

d-Dates; e-Inclination; f -Authorities; 
1-Hornsound; 2-Ice Sea; 3-Swarte-Clif f 
(New Zero); 4-Cape-North; 5-Coast of 
Norwary; 6-Baffin Bay; 7-Coast of Norway; 
8-Frobisher Strait; 9-Near Shetland; 
10-Uranibourg; 11-London; 12-Rouen; 
13-Paris; 14-Dauff beuren; 15-Tours; 
16-Ferrare; 17 -Rome; 18-Naples; 19-Lisbon; 
20-Messina; 21-Palerma; 22-Malta. 

23-June; 24-May; 25-July; 26-November. 
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For the inclination in the XVI th and XVII th > 
consult the following publications: 


-Hansteen, Untersuchungen liber den Magnetismus der Erde , 
t . I , p . 38 . 

-For tun io Affaitati, De peculiari magnetis ad polum 
descensu etc . (Physicae et astronomicae considerationes, 
Venet, 1549.).. 

-R. Norman, The New Attractive , London 1580. 

-Cabaeus, Philo sphia magnetica , Ferraiae 1629. 

-Kircher, Magnes, sive de Arte magnetica , p. 401, 424 
(.1641 edition) . 

-Jacobus Grandamicus, Noua demonstratio Immobiljatis 
Terroe , cap. 4, sect.l, La FlSche 1645. 


Tycho-Brahe' s and Kepler's observations were made 
from the position of the point of convergence of the 
northern light beams observed by these astronomers. 

Tycho-Brahe 1 s meteorological log-book recorded 
at Uranibourg from 1582-1597, published in Copenhagen in 
1876 by Paul la Cour, p. XLV. 

-Marian, Treatise On The Northern Lights , p. 130, 188 
(Continuation of Memoires de l'Academie Royale des 
Sciences, 1731) . 


1 2 . 


The objective of our work, was primarily to offer a 
numerical theory of the earth' s magnetic attraction, ap- 
plicable to any time, and we have arrived at a solution 
to this problem by approximating the truth. 

This solution, although it already appeared in an 
excessively simply analytical form, was nevertheless of 
a purely empirical nature. And this type of 'solution is 
not satisfactory to the scientific spirit. The true ob- 
jective of the theory is to trace back to the physical 
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causes and to thus give the formulas their theoretical 
justification, by starting with the principles of physics 
or mechanics alone. 

From this standpoint, the approximated solution found 
by inductive reasoning, may be considered an approximatively 
known integral of the problem. To arrive at this objective, 
it is necessary to deduce from it the forces exerted at 
each moment. This will make it possible to express the 
differential equation which the question is dependent upon, 
for ithe hypothesis adopted. 

In the present chapter, we will take up this difficult 
problem, by slightly exceeding the boundaries of precision 
which v;e scrupulously tried to observe until now. 

Let us state first that any complete theory of the 
earth's magnetic attraction seems possible only if the 
effects observed come from a constant cause which always 
acts in the same way. If, for example, the changes were 
due to a slow cooling of the earth, or to the transfer of 
matter to unknown depths, how can they be as constant as 
they are? 

Let us therefore assume that the changes are due to 
a constant and always uniform cause, and above all to a 
continuous reraagnetization of the earth, which is assumed 
to be a continous magnet. Now we will ask ourselves: 
what could a magnetic field be, at each'dnstant, which 
would produce the remagnetization observed? 

Let us assume there is a sphere magnetized by any, 
but variable, mode. Let V be the sphere’s potential for 
an outside point at any moment. Also let ft be the poten- 
tial at the same point due to the induced magnetization 
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developed within the sphere, per unit of time, 
by a system of electrical currents or by any other for- 
eign cause. By borrowing a mode of language specific to 
mechanics, we will call this the virtual potential of the 
sphere. Once this is established, the total or instantan- /3Q 
eous potential V will be related to the virtual potential 
by a relationship expressed: 



« ft 


where ft is estimated to contain a numerical factor which 
is equal to the quantity of the continuous magnetism dev 
eloped within the sphere per unit of time. 


We are proposing to determine this virtual magneti- 
zation as a function of the instantaneous magnetization, 
so as to satisfy the conditions created by the earth. 

We will consider an especially simple case, that of 
a homogeneous spherical layer, or of a sphere composed of 
concentric and homgeneous layers. 


The potential ft of the magnetization induced in the 

sphere is assumed to be caused by an outside magnetic 

field V' which we are trying to determine. Poisson's 

theory told us to determine ft as a function of V. on the 

assumption that ft is expanded into a series of spherical 

(i) 

functions, where each ft ‘ term of order i is directly 

(i) 

proportional to a term of the same V' order in the 
inductor field V r , so that we may express: 


by letting N. be the factor of proportionality. 

'I* ^ 
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We are therefore assuming that the continuous mag- 
netism follows the same laws as a temporary magnetism 
aroused by induction. The time factor may be estimated 
to be contained in the factor N.j which also contains 

the outside and inside half-diameters a, b of the hollow 

(i) 

sphere, the order i of the function V » and the coef- 
ficient of magnetization induced k in the middle. 

An analogous expression 'is obviously appropriate for 
the case of a sphere composed of concentric and homogen- 
eous layers, or even for a sphere where the coefficient 
of magnetization is a function of the vector radius alone. 
We simply have to replace the expression above by a sum, 
where each element corresponds to an isolated layer, or 
by an integral. 

By letting N 1 . be the factor of proportionality and 

* l * 

by setting: 


£2 t £2^°^ + £2 ^ ^ + £2^^ + etc., 
we will always have the following relationship: 

£2^ = N 

Once this is established, the differential equation 
which defines the total potential of the earth V at each 
instant, as a function of time, will be: 

= N V ,fo) + N 1 V ,(1) + N 8 V ,(2) + eto. 3 

N o jN 1 jN 2* var y from one term to another 


dV_ 

dt 

where the numbers 
in the series. 
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If the second member is zero, this equation becomes: 


V « Arbitrary function of y , tar; 


the arbitrary function being determined by the initial 
sta te. 

To solve the proposed equation in the case where the 
second member is not zero, we expand V into a series of 
Laplace functions and we will apply the method of varia- 
tion of arbitrary constants to it. We will therefore set: 

V = F C0) + F<‘> + F w + etc.; V (l) = Fj> + V? + ...+ Vf ; 

K = ' • «? cos « ( « + 


and the analogous expressions for the V' > which are 

Yl / * i f * \ 

deduced from the first ones by replacing a' ^ i and 3' ^ 

with a (%) and 

n n 

(i) (i) 

If we vary the s and 6 n s, and identify the 
cosnw and sinmzr coefficients of the two members,, we will 
have a system of conditions [we have omitted the common 
factor : (/t — /**)"—] •' 


C0S rf -fif - na n sin n Pn = Xt «? cos nff , 
r/a f0 ,7/9(0 

siu + Wft (;> cos nff = ^ sin ; 


from which we may promptly deduce: 


*«> j 4| -y '(0 

f - ^ CM W = 'W? ~ tf) 
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poor qualit y 


If a 


(i) 


¥1 /• \ 

in the observations, it is necessary for n ( 3 w 

w 


does not contain a secular term, as shown 


& (i> ) 

p n ' 


(2k + 1)~, by letting k be any whole number. On this 
assumption, the second equation is expressed: 




(i) 


n_ 

dt 


= + 


N . a 
t n 


,(i) 


na 


(i) 

n 


If it is integrated, it becomes: 


= ± N.a’ (i) 

n t n . , ( i ) 

(i) + Y n - 


n a 


(i) 


n 


and by letting y^ be an integration constant. 


The observations show that the factor t is positive, 
and since is negative, it is therefore necessary to 
use the lower sign, which shows that k is an odd number. 
It follows that for any number 

71 71 a 

of full circumferences. 


The conditions established by the data will there- /3 2 

fore be satisfied. 

The variations observed of the earth's magnetization 
are therefore caused by a virtual magnetic field related 
to a primitive field by the following conditions: 


"The modules of the virtual field are very small 
fractions of the modules of the instantaneous field, where 
the arguments differ from these by the angle 5- . " 


On this assumption, the average yearly motion m 

N a d ) 

of any given term is equal to: — _i n and the 


na 


(i) 

n 


(i) 

n 
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corresponding magnetic moment's ■piiye of rotation about 
the earth's axis will be 2 ’ ]rna n 

i - n 

If we let be the yearly motion of we may 

set: 




t + y 



( * 7 , ) 

Y is the initial longitude which corresponds to the 
1 1 

origin of time. 


13 . 


Following these preliminary considerations, we will 
put the question in specific terms and we will ask our- 
selves which forces are capable of exerting the required 
virtual magnetization, for each instant? 

It is not our intention to dispute the countless 
hypotheses put forth in this regard. We simply wish to 
propose What we consider to be the most likely hypothesis. 

In our opinion, the phenomenon of the secular varia- 
tions of the earth's magnetization is a large phenomenon 
of electromagnetic induction taking place in the rarified 
layers of the atmosphere, set in motion by the earth's ro- 
tation.' Below, u we will see that this hypothesis effectively 
accounts for the essential traits of the phenomenon. 

The existence of induction currents in the atmos- 
phere seems well established by other circumstances. 

The role of the terrestrial atmosphere in the diurnal 
oscillations of the magnetic elements has been well-known 
by scientists since Schuster's report on this subject 
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which he had inserted in volume CLXXX of the Philosphical 

Transactions in London. The result is that this variation 

(o ) 

corresponds to a term of form Y^ of the force function, 
where the argument increases uniformly with time and de- 
scribes a full circumference in 24 hours. This term seems 
to be sort of related to the diurnal oscillations of the 
barometer. 

X have demonstrated myself, in an unpublished note, 
that the semi-diurnal lunar variation of the terrestrial 
magnetic field is caused by electrical currents circula- 
ting in high regions of the atmosphere. The corresponding 
term in V is the expression Y^ J I demonstrated next 
that an atmospheric tidal wave would generate a term of 
this form. 

In the cases mentioned, it is due to periodic motions /3 3 
that the atmosphere becomes the center of induction cur- 
rections under the action of the terrestrial magnet. More- 
over, information is not lacking that there is :.a uniform and ro- 
tary motion of the high layers of the atmosphere, which 
seems to play an important role in the phenomenon of polar 
light . 

Most likely, there is nothing a priori in this sup- 
position that the rarified layers of the atmosphere become 
the center of electrical currents, capable of reacting to 
the earth itself, and which thus produce this continuous 
remagnetization of the earth which we observe the effects of. 

14. 


This brings us to study the currents induced in a 
homogeneous spherical layer whose particles are animated 
in any given rotary or oscillating motions. 
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First, it is easy to see that the oscillating motions 
having a speed potential will generate a system of induc- 
tion currents in such a manner that the current function 
will be of a higher order than that of the term that 
generates it and it will never have the same expression, 
as required by the observations. 

This is no longer the case for the rotary motions 
of the atmosphere. Consequently, we will not discuss the 
problem in all its generalities, butwill examine the hypo- 
thesis which seems to fit the facts the best, and this is 
the case of a spherical layer undergoing a uniform rota- 
tion about a stationary axis. 

Moreover, this problem was the subject matter of 
Her z 1 first work ( OEuvres de Herz , t. I), so that we may 
benefit from the results of his analysis. 

Let us consider a thin spherical layer, of radius i ? ' 
undergoing a uniform rotary motion, of angular velocity 
to, about the axis of a concentric spherical magnet. Each 
term 



of the function of the magnetic forces, generates a sim- 
ilar term in the running function cj) which has the follow- 
ing expression: 


B co 
i v. 



sin (>M + H/jjl) , 


where x represents the specific resistance of the surface.* 


*Remember that we have not defined the magnetic potential 
V = — / contrary to most modern authors who have 
adopted the notation V - f~. 
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This system of induced currents itself generates 
a magnetic field, within the layer, and under our hypo- 
thesis, this magnetic field is identical to that V ' 
discussed above. 


Let us evaluate this field. By borrowing Maxwell’s 
formula in his Treatise (volume II, art. 67 2)., we find, 
for an inside point of distance v from the center of the 
sphere, L+l. !£ £(')' (f,f .Pj'. ,(»» + .<")• 

This is indeed the expression of the required form. 
The hypothesis is therefore completely verified. 
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By expressing the expression for the inductor field 
in its simplest form: 

P? cos (nt7 -f n/"*) , 


a ) 

we find, for argument 0'^ *, the following relationship 

(i) /VI TT 

n.Y n “ n $ n ~ ~ 2 * If w is negative? this will bring 

us to: cos (jittr + wjj'lf) = + sin (wnr + and we must therefore 

take the lower sign and set: 


o __ j <> 77 

rn ? i 2H 


Accordingly, if the relative motion of the upper 
layers of the atmosphere is directed in the opposite dir- 
ection to the diurnal motion, the motions observed of 
a \ 

0^ coincide with the theory. 


(i) 

For the module , we have the following expres- 


sion: 


*os is positively accounted for in the direction of 
increasing angular velocities. 
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Above, we found that (no, 12) the yearly motion of 3 


(i) 


is equal to n 


tCii 


n 


n a 


Til 


; for functions of different orders, 


n 


while the yearly motion is dependent upon both the order 
i and N^, AL being a complication function of i and on 
the law of variation of the induction coefficient within 
the earth. In >our case, the factor is essentially in- 
dependent from ij the yearly motion should be proportional, 

for the terms of different orders, to i+1 , and there- 

(U+l)i ' 

fore, specifically, for the functions of the first four 

orders, it would be approximately within the ratio of the 

fractions 2 , 3,4 5 . This condition is ap- 

371 5.2 773 ' 974 

proximately fulfilled, with one importnat exception, by 
the motions which the observations enable us to establish 
with some certainty. 


15. 


The preceding analysis accounts for the main part 
of the changes in the earth's magnetization, and nothing 
proves that these are exact laws of the phenomenon. Con- 
versely, it is probable that the changes are actually more 
complex, and it would be useful to become aware of what 
we might expect. 


The formulas which we have given are appropriate for 

the case of a homogeneous sphere or a sphere composed of 

concentric or homogeneous layers. A heterogeneity of the 

(i) (i) 

magnetized sphere introduces periodical terms in a and 8 
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in the solution. 


Let us quickly examine what happens in a heterogen- 
eous sphere which deviates very little from homogeneity. 

( d* ) 

As a first approximation, we may still assume the a s 
are constant and the are in the form 

Yl 7> 


n 


Ci) 


The differential equation which each Y n term of 


the function of forces satisfies, will still be: 

dr!? 


~w 


3®. 


a ) 

where Z„ is expressed: 


n 


cos » (tar + (j'f) . 


However, here, a* and ^ will contain the periodical 

, . n n 

terms in 2 n&Y , or, which amounts to the same thing, m 

( * ) c dsf a ) 

2nm { ^ t 3 m n being the yearly motion of ^ • 

(%) (i) cos wo’ 

By expanding Y ^ and Z^ and identifying the 

and sinm * coefficients of the two .'members , we will again 

have, by virtue of a known theorem: 


d&° 


W cos n *»? “n sin cos n/f® , 


da!? 


ir sinM &° + COS = Jr, a<° sin nfl?; 


from which we may easily deduce: 

'^=y,J®co S »(«'-©; -a?- 


net 


m sin w (?? - $?) 




:«) 


In these equations and /S' ( n ° , and consequently 

and jjF? have received small periodical increases da ( n °, d/jJ?, 
etc. which we are trying to determine. 
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By deducing the equations, which are analogous for 
the case studied earlier, from these equations, we will 
have the equations with the following differences: 


I dd«P id K,) a' m \ 


cUff" N. I a w a' {,) \ 

i rr “ * (> + J ■#>) cos * OW - « - 


<(/) 


* r ( < 


jo* 




If the sphere under consideration does not deviate 
much from homogeneity the angle will be close 

will differ very little from 


to 


re 


m 


and the ratio 


;-«) 

Li 

rw 




. By expanding the second members according to 
the powers of very small quantities ? and t defined by the 


it 


and 


;'«>■ 
L n 

;«) : 








\ (i + r), and 


relationships o = w — ■ p\‘) + 

by omitting the fourth order quantities, we will have, 
by inverting the order of the differentiations in the 
first members: 
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'^r + «-?> 


M ■ 


£* 5; r ; 

n „«> ' 




Without going to the bottom of the question for 
the time being, we can predict that the quantities a 


and t will be periodical with respect to 2w3 


U) 


n 


for 


reasons of symmetry, and consequently with respect to 
C *1* ) 

2nm „ 't. We therefore have: 
n 


da 


(il 


n 


dt 


d 3 

and 6 dt 


(i) 


til 

= (Periodic terms in 2nm t. 


n 


til 


These equations will cause the a s and 3'' w s to become 

n n 

integrated with, the trigonometric series proceding accord- 

(i) 

i.ng to the multiples of the angle 2nm t. The periodic 


n 


time of the first term of the series will be 


il . , 

2nm t ' 1 ' 
n 
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a ,( '> 

j a 

We can easily see that by assuming «w and » 
are of the same infinitesimal order, 
will be of the order of V„°> This will occur, for 
example, if the heterogeneous sphere is assimilated with 
a homogeneous sphere surrounded by a circular ring pass- 
ing bhrough its poles. Thus, the periodic terms in « 
and $ will be of the same order. 

By closing this superficial discussion here, we 
may summarize the results as follows: 

"In a heterogeneous sphere which substantially ap- 
proaches a homogeneous sphere, the heterogeneity intro- 

(i) 

duces periodical terms into the expressions of the 

(i, ) 

and 3^ whose cycles will be the half-periodical times 
of the corresponding term in the function of forces, and 
its multiples" . 

These are the main special characteristics which 
we might hope to see come out of a rriore thorough discus- 
sion of the overall observations made, if the ideas put 
forth in this report are accurate. 


> , • 
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